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relative  to  that  observed  when  pH  3.0  HNO3  was  applied.  All 
other  treatments  showed  net  losses  of  cations.  Net  losses 
occurred  for  the  pH  3.7  treatments.  Apparent  net  losses 
were  greatest  for  the  pH  4.7  treatments,  small  net  losses 
from  the  topsoil  leached  with  pH  3.0  H2SO4  were  likely  due 
to  overestimation  of  the  charge  borne  by  Al.  Retention  or 
loss  of  cations  appeared  related  to  SOj*'  sorption- 
desorption. 

A simple  mixing-cell  transport  model  coupled  with  a 
chemical  equilibrium  submodel  that  included  anion  sorption 
was  used  to  describe  Rt-induced  cation  leaching.  Although 
the  model  predicted  later  cation  breakthroughs  for  leaching 
with  H28O4  than  with  HNO3,  predicted  breakthroughs  were 
later  than  experimentally  observed.  Due  to  wide  variations 
in  soil  pH,  the  Freundllch  isotherms  used  were  inadequate  to 
describe  504^*  sorption.  Modification  of  the  chemical 
submodel  to  include  time-dependent  cation  exchange  and 
explicit  pH-dependence  of  anion  sorption  are  apparently 
needed  to  achieve  agreement  of  model  results  with  experi- 
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IKTRODUCTIOH 


Soil  acidification  reaults  Cron  any  regime  of  procesees 
that  results  in  net  removal  of  basic  and  addition  of  acidic 
cations.  The  natural  process  of  leaching  that  tends  to 
occur  as  water  moves  through  soil  nay  effect  soil  acidifi- 
cation if  the  losses  of  bases  are  hot  offset  by  inputs  by 
atmospheric  deposition  or  nutrient  recycling  or  releases  by 
mineral  weathering.  Biologic  processes  within  the  soil  tend 
to  generate  H'*’  ions,  oxidation  of  H and  S liberates  and 
respiration  releases  COj  which  combines  with  H2O  to  form  the 
weak  acid.  H2CO3.  Some  nutrient  uptake  by  plants  is  also 


acidifying. 

Natural  acidifying  processes  may  be  aggravated  by  soil 
management  practices.  Nutrient  recycling  is  interrupted  by 
removal  of  ions  during  harvesting.  Amendments  of  oxidisable 
K fertilizers  tend  to  accelerate  H'*'  production  in  the  soil. 

Soil  acidification  is  also  effected  by  anthropogenic 
(i.e.,  man-made)  deposition  of  sulfuric  and  nitric  acids. 

The  effect  of  this  additional  acidifying  process,  however, 
has  been  estimated  to  be  small  compared  to  that  due  to 
amendments  of  N fertilizers  (HcFee  et  al.,  1$83).  Further- 
more. management  practioes  of  agricultural  soils  regularly 
include  application  of  limestone.  HcFee  et  al.  (1983) 
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concluded  that  for  an 
effect  of  acid  dspositicn  on  base  leaching  and  acidiflca*- 
tion,  particularly  the  aobilizaticn  of  Al,  may  be  dele- 
terious to  plant  grovth. 

One  factor  which  would  seen  to  nltlgate  these  poten- 
tially adverse  effects  is  sorption  of  anions  from  the  acidic 
leaching  solution.  If  input  anion  species  are  sorbed  to  a 
greater  extent  than  native  anions  are  desorbed,  cation 
leaching  must  be  retarded.  Furthermore,  if  different  anions 
exhibit  different  tendencies  for  sorption,  the  anionic 
composition  of  the  input  solution  would  affect  cation 
leaching  rates.  Thomas  (1960)  has  demonstrated  such  effects 
for  anion-sorbing  soils. 

The  anion  species,  MO3”  and  dominant  compo- 

nents in  acid  rain  [Bubenick  et  al.,  1983;  Jac)cs  et  al., 

1984) . Differences  in  the  extent  of  sorption  of  these 
species  would  effect  differences  in  cation  leaching,  acidi- 
fication, and  mobilization  of  Al  caused  by  h*  loading.  The 
principal  objective  of  this  investigation  was  to  discern 
differences  in  leaching  cates  for  several  common  cation 
species  from  an  anion-sorbing  soil  leached  with  either  HNO3 
or  H2SO4.  Since  pH  is  known  to  affect  anion  sorption,  these 
acids  were  applied  at  three  different  pHs.  Another  objec- 
tive was  to  evaluate  a predictive  model  for  describing  the 
effects  of  pH  and  anion  sorption  upon  cation  leaching.  The 
intent  here  was  to  assess  the  applicability  of  simple 
transport  and  chemical  equilibrium  models  for  the 


conplicated  system  of  physical  and  chemical  processes  which 
occur  during  Che  transport  of  acid  solutions  through  soil. 
The  transport  model  used  was  a nlxing-cell  model.  This 
model  simulates  hydrodynamic  dispersion  by  a simple 
algorithm.  Accordingly,  computational  time  is  reduced  and 
the  code  is  simple.  The  principal  advantage  of  the  chemical 
equilibrium  model  used  is  that  cation  exchange  selectivities 
are  estimated  from  solution  and  sorbed  phase  concentrations 
of  the  unaltered  soil.  Detailed  studies  of  binary  exchange 
reactions  are  not  required. 
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LITERATURE  REVIEW 

Thonas  (I960)  compared  ewperinencally  determined 
breakthrough  curves  (BTCs)  for  concentrations  of  ca2+  In 
effluent  from  columns  of  Ca2*-saturated  Volo,  Cecil,  Lloyd, 
and  Davidson  series  soils,  leached  with  various  Al^'*'  and  K* 
salts,  with  predicted  BTCs  based  on  a aimple  chromatographic 
theory.  Agreement  was  good  for  the  yolo  soil  but  poor  for 
the  three  Piedmont  soils  (Cecil,  Lloyd,  and  Davidson).  For 
the  latter  group  of  soils,  actual  breakthrough  for  Ca^'*'  was 
later  than  predicted.  These  deviations  were  attributed  to 
electrolytic  imbibition  which  is  the  combined  retention  of 
anions  and  cations.  That  author  also  investigated  anionic 
effects  on  electrolytic  imbibition.  While  Cl“  and  NOy" 
salts  effected  some  retardation  of  Ca^'*'  with  respect  to 
theory,  SOa^”  salts  resulted  in  greater  retardation  for  ca^'*' 
movement.  Retardation  was  greatest  when  H2PO4'  salts  were 
used.  Chloride  and  NO3”  anions  tend  to  be  sorbed  by  soil  to 
a lesser  extent  than  are  SO42*  or  H2PO,-  anions.  The 
results  of  this  study  have  en  important  implication  for  the 
problem  of  accelerated  leaching  of  cations  from  soils 
subjected  to  acid  deposition.  For  soils  that  exhibit  anion 
sorption,  the  nature  and  relative  amounts  of  the  anions 
deposited  will  affect  the  rate  of  H+-induced  cation 
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leaching.  According  to  Johnson  et  si 
can  be  acidified  by  infiltrating  acid  rain  without  being 

infflobllited  in  the  soil.  Lee  (198S)  irrigated  lysineter 
boxes  containing  Bohannon  (Typic  Halpunbrept)  soil,  a SOa^" 
sorbing  soil,  and  planted  with  either  young  red  alder  (Alnua 
rubra  Bong.)  or  sugar  maple  (Acer  saccharun  Harsh.)  with 
H2SO4  solutions  of  several  pHs.  After  three-and-one-half 
years,  Ca^***  and  Hg^*^  had  been  leaohed  to  a depth  of  only  15 
to  30  cm  for  the  pH  3.0  treatment.  This  depth  coincided 
with  that  at  which  S04^~  was  being  removed  from  solution  by 
sorption.  Huete  and  HcColl  (1984)  found  that  smaller 
amounts'of  cations  were  leached  from  an  Alfisol  subsoil  and 
an  Oxisol  topsoil  when  H2SO4,  rather  than  HNOj,  was  applied. 

Cation  mobility  is  clearly  dependent  upon  anion 
mobility.  Accordingly,  this  study  was  underta)cen  to  deter- 
mine the  effects  of  differences  in  the  mobilities  of  NO3* 
and  SOa^”,  components  in  acid  rain,  upon  cation  leaching 
rates  from  an  anionic  sorbing  soil.  Factors  affecting  the 
mobilities  of  these  anions  are  subsequently  discussed. 

Factors  Affecting  Anion  Mobility 
Anion  movement  is  determined  by  the  combined  effects  of 
convective-dispersive  transport  and  sorption-desorption 
processes.  For  1-dinensional  transport,  reversible  sorp- 
tion, and  steady  water  flow  the  transport  equation  for  a 
given  species  may  be  expressed  as 


where  C is  solute  concentrstion  (mole  S is  sorbed 

phase  concentration  (mole  P is  a sink  or  source  term 

(mole  m*3  a"lj  for  removal  or  addition  of  a species  by 
precipitation/dissolution,  irreversible  sorption,  or  micro- 
bial activity,  9 is  volumetric  water  content,  o is  bulk 
density  (kg  m"3),  D is  the  hydrodynamic  dispersion  coeffi- 
cient (m^  s“i),  V is  the  pore  water  velocity  (m  s**^) , X is 
distance  (m) , and  t is  time  (s). 

Generally  the  greater  the  extent  of  sorption,  the 
greater  is  the  retardation  of  an  anionic  sorbing  species 
relative  to  a non-sorbing  species.  Several  factors  influ- 
ence the  amount  sorbed  per  unit  mass  of  soil  as  well  as  its 
time  rate  of  change.  Major  factors  affecting  SO^^-  sorption 
include  1)  solution  concentration,  2}  pH,  3)  nature  of 
sorbent,  and  4)  nature  of  other  ionic  species  present 
[Gaston  et  al.,  1996).  These  factors  likewise  influence 
NOj"  sorption  and,  consequently,  its  mobility. 

sulfate  sorption  has  been  observed  variously  to  either 
conform  to  the  Langmuir  isotherm. 


[2] 


In  which  a maxlnum  In  sorption  is  reached  with  increasing 
solution  concentration,  or  to  the  Freundlioh  isothem, 


C3] 


in  which  no  sorption  maxinum  occurs.  In  these  equations  C 
is  solute  concentration  (nole  n~^),  S is  sorbed  phase 
concentration  (mole  kg“i) , b is  the  sorption  maxinun  (nole 
kg“i)  and  K Bole“i) , k (inolel”^  kg"lj  and  n are 
empirical  constants.  Kinjo  and  Pratt  (1971)  found  the 
Langmuir  equation,  (2],  better  described  HOj"  sorption  for 
the  subsoils  of  two  Andepts,  an  Alfisol,  and  an  Oxisol. 

Models  used  for  predicting  effects  of  acid  deposition 
on  soil  chemistry  have  commonly  assumed  reversible  isotherms 
for  SO42*  sorption,  sulfate  sorption  may  not  always  be 
reversible,  however  (Aylmore  et  al.,  19£6] . Apparent  non- 
singularity due  to  irreversible  sorption  in  SOa^"  isotherms 
would  affect  leaching  rates  for  soil  cations  in  the  situa- 
tion where  heavy  inputs  of  SO42-  to  a soil  were  followed  by 
low  levels  of  SO42-  input.  As  pointed  out  by  Reuss  and 
Johnson  (19S6) , non-singularity  in  this  case  would  diminish 
H7‘-induced  leaching  of  cations. 

Sulfate  sorption  generally  increases  with  decreasing 
pH.  nitrate  sorption  behaves  similarly  (Kinjo  i Pratt, 

1971)  Black  s Haring,  1979).  Although  this  behavior  is 
consistent  with  the  increase  in  positive  sites  for  anion 
adsorption  with  decrsasing  pH  and  accounts  for  the 


effect 


pH  on  NOj"  sorption  (Hingston  et  si.,  1973),  304^“  adsorp- 
tion say  also  occur  at  nautral  aitas  (Rajan,  1978;  1979). 

In  thsss  invsstigations,  Rajan  observed  an  increase  in 
surface  negative  charge  due  to  304^*  adsorption.  Couto  et 
al.  (1979)  also  found  that  SO**”  sorption  slightly  increased 
the  CSC.  Kingston  et  al.  (1973),  however,  clais  that  no 
evidence  exists  for  SO**”  adsorption  at  neutral  sites, 
sulfate  adsorption  proceeded  with  neutralization  of  positive 
sites  so  that  the  Unit  of  SOa*”  adsorption  at  a given  pH 
of  non-specifically 
■ sorption  occurs  only  by 


adsorption  and,  on  an  equivalent  basis,  cannot  exceed  the 
naxinujD  of  non-specifically  adsorbed  anions,  then  cations 
would  be  leached  fron  an  anion-sorbing  soil  just  as  effec- 
tively by  H2SO4  as  by  HNO3.  As  Singh  et  al.  (1980)  pointed 
out,  imnobillzation  of  input  anion  will  result  in  retarded 
leaching  of  cations  only  if  a less  than  equivalent  amount  of 
native  soil  anions  is  simultaneously  displaced.  In  light  of 
the  leaching  studies  by  Thomas  (I960)  and  Huete  and  HcColl 
(1984),  which  demonstrated  retarded  cation  losses  when 
SO**”,  rather  than  the  non-spacificslly  adsorbed  HO3”,  was 
the  major  anion  in  the  system,  the  mechanism  proposed  by 
Kingston  et  al,  (1972)  is  not  solely  responsible  for  SO4*” 
sorption.  (Proposed  mechanisms  for  SO4*”  sorption  are  dis- 
cussed later.) 

Among  inorganic  sorbents,  aluminum  and  iron  oxide 
materials  and  kaolinitic  minerals  exhibit  the  greatest 


capacity  to  sorb  S04^~  (Gaston  at  al.,  1986).  These  same 
inorganic  materials  were  responsible  for  sorption 

(Klnjo  i Pratt,  1971:  Black  S Waring,  1979).  The  effect  of 
organio  matter  on  S04^~  sorption  is  unclear.  Although 
several  authors  found  that  organic  matter  contributed  to 
SO^^*  sorption,  some  found  that  it  exerted  a negative 
influence  on  sorption  (Gaston  et  al.,  1986).  Black 

and  Haring  (1979)  found  that  organic  natter  tended  to  reduce 
N03~  sorption.  According  to  these  studies,  anionic  effects 
on  cation  leaching  should  be  expected  for  soils  containing 
appreciable  amounts  of  1:1  clays  (kaolinite,  for  example) 
and  aluminum  and  iron  oxide  materials. 

Competition  of  other  sorbing  anion  species  for  sorption 
sites  tends  to  reduce  the  extent  of  SO^^”  sorption.  Phos- 
phate has  been  shown  to  decrease  so^*-  sorption  (Ensmlnger, 
1954;  Kanprath  et  al.,  1956:  Chao,  1964).  Chloride  (Chang  6 
Thomas.  1963;  Cebbardt  s Coleman,  1974)  and  NO]*  species 
(Chao,  1964)  are  poorly  competitive  with  3042*.  Khanna  and 
Beese  (1978)  in  Hast  Germany  shoved  that  a pulse  application 
of  cl*  and  SO]*  to  columns  of  an  acid  podaolic  brown  soil 
saturated  with  SO42-  did  not  desorb  3042*.  It  may  be 
conversely  concluded  that  SO]”  adsorption  is  restricted  in 
presence  of  3042-,  Preferential  sorption  of  SO42*  is  at 
least  partially  understood  in  terms  of  differences  in 
surface  binding  of  these  species.  Nitrate  is  nonspecifi- 
caily  bound  to  positive  surface  sites  (Hingston  et  al., 

1972)  by  Coulombio  attraction,  sulfate,  an  the  other  hand. 


Is  specifically  bound  to  surface  U or  Fe  by  displacenent  of 
HjO  (Kingston  et  al..  1972)  and/or  OH“  (Bajan,  1978  and 
1979).  Serna  et  al.  (1977)  presented  spectroscopic  evidence 
for  such  weak  binding  for  NO3*  compared  to  S042*. 

In  addition  to  anionic  affects  on  so^*-  sorption,  Chao 
et  al.  (1963)  demonstrated  that  the  composition  of  the 
cationic  milieu  affects  S04^~  sorption.  The  higher  the 
valence  of  the  predominant  cation  in  the  system,  the  greater 
was  the  extent  on  S04^~  sorption.  The  authors  suggested 
this  effect  could  be  due  to  either  the  formatlcn  of  a 
bridging  complex  in  which  a trivalent  or  divalent  cation 
links  304^*  to  a negative  exchange  site  or  to  reduced  anion 
repulsion  when  the  diffuse  double  layer  contains  such 
cations. 

In  addition  to  the  four  major  factors  discussed  above, 
S04^~  sorption  ie  also  affected  by  time  and  temperature. 
Rajan  (197S)  found  that  although  S0^^~  sorption  onto  hydrous 
alumina  was  about  $0%  complete  within  10  min.  reaction  time, 
equilibriiun  was  not  reached  until  three  hours.  Chang  and 
Thomas  (1963)  found  that  soa^*  sorption  on  2 Cecil  soils  and 
a Georgia  kaolinite  coated  with  Fe  oxide  continued  over 
several  days.  Barrow  and  Shaw  (1977)  found  that  the  removal 
of  3042"  from  solution,  availability  of  SOa*"  to  plants,  and 
desorption  each  decreased  with  increasing  contact  time 
between  solution  and  sorbed  phases.  Singh  (19B4)  described 
the  kinetics  of  3042“  sorption  for  B horizon  material  of  a 


Typlc  Udlpsamiitent  and  a Typic  DystrochrepC  using  a 2-rata 
ccmscann  and  multiple-order  rate  equations. 

On  the  basis  of  thermodynamics,  anion  sorption  should 
be  affected  by  changes  in  temperature.  Barrov  and  Shaw 
(1977)  and  Singh  (19S4)  shewed  that  sorption  increased 

slightly  with  increasing  temperatures  near  ambient  tempera- 


on  the  basis  of  these  factors  affecting  NO^'  and 
S04^"  sorption  it  can  be  generally  concluded  that  the 
mobility  of  these  anions  will  be  least  (i.e.,  sorption 
should  be  greatest)  in  a soil  composed  of  1:1  clays  or 
aluminum  and  iron  oxide  materials,  having  a distinctly 
acid  pH,  and  containing  low  amounts  of  a preferentially 
sorbed  anion, 

Results  of  studies  by  Thomas  (I960)  and  Huete  and 
Mccoll  (1984)  Indicate  that  for  soils  in  which  a difference 
in  mobilities  of  these  anions  exists  NO3*  is  more  mobile 


than  SO42-.  The  lower  mobility  of  occurs  because,  for 

the  same  sorbent  materiel  and  pH,  a greater  portion  of  the 
total  amount  of  the  anion  exists  in  the  sorbed  phase.  The 
machanisBs  whereby  S0,2-  is  retained  by  the  soil,  compared 
to  that  by  which  HOs*  is  adsorbed,  account  for  this  greater 
sorption  of  sOa^”. 


Sorption  Mechanisms 

Sorption  of  NO3*  is  by  Coulombic  attraction  to  posi- 
tively charged  sites  on  aluminum  and  iron  oxides  and 


Icaolinitic  minerals  (Kingston  et  al.,  1972).  Sorption  of 
3042"  may  occur  by  any  of  three  general  mechanisras  (Harvard 
4 Relsenauri  1966) . These  are  1)  coordination  with  aluminum 
or  Iron  in  oxide  materials  via  displacement  of  coordinated 
H2O  or  0H~  and  a elmllar  exchange  on  the  edges  of  1:1  type 
clays,  2)  a mechanism  whereby  304^*  and  cation  are  both 
retained  (termed  "molecular  adsorption"),  and  3)  chemical 
precipitation  reactions. 

In  the  ligand  exchange  mechanism,  S04^~  may  be  adsorbed 
as  a monodentate  ligand  displacing  HgO  from  positive  sites 
(Kingston  et  al.,  1967;  1972).  Sulfate  nay  also  be  adsorbed 
as  binuclear  bridging  complex  displacing  H2O  or  OH~  from 
adjacent  metal  atoms  (Bajan.  1978;  1979).  Parfitt  and  Smart 
(1978)  presented  spectroscopic  evidence  for  bidentate  804^* 
adsorption.  Their  overall  reaction  for  SO4*'  adsorption  on 
iron  oxide  materials,  analogous  to  that  proposed  by  Rajan, 
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Molecular  adsorption  may  be  conceptually  viewed  as 
concurrent  specific  adsorption  and  cation  exchange  (Harvard 
4 Relsenaur,  1966).  According  to  Chang  and  Thomas  (1963), 
an  adsorbed  acidic  cation  hydrolyzes,  releasing  H"*,  allowing 
for  the  adsorption  of  SO42-  co-ion,  and  SO42-  displaces  OH". 
Neutralization  of  H"'  and  OH"  drives  the  reaction.  Another 


form  of  molecular  adsorption  is  the  physical  entrapment  of 
S04^~  salts  in  soil  micropores  (Harvard  and  Reisenaur, 

Several  authors  have  suggested  that,  while  is 

adsorbed  by  a host  of  materials,  it  may  also  be  precipitated 
in  the  presence  of  these  materials  as  basaluminite, 
Al4(OH)ioS04'SH20,  or  alunite,  KAI3 (OH) g (SO4) 2 (Caston  et 
al.,  19se).  Whether  these  minerals  control  soil  solution 
504^~  is  unclear.  Adams  and  Rawajfih  (1977)  found  that 
basaluminite  might  be  controlling  S04^~  solubility  in  a 
K2SO4  amended  sandy  loam.  The  data  of  Holt  (1981)  suggest 
that  alunite,  in  conjunction  with  gibbsite  or  kaolinite, 
might  control  3042*  in  solution.  Khanna  and  Basse  (1978) 
could  not  demonstrate  that  either  mineral  was  controlling 
the  solution  concentration  of  304^*.  Gebbardt  and  Coleman 
(1974)  found  that  in  the  soil  solution  of  B horizon  material 
of  an  Andept  ion  products  for  basaluminite  spanned  a range 
from  supersaturation  to  undersaturation. 

These  mechanisms  by  which  S04^~  is  sorbed  render  it 
more  tightly  and  preferentially  bound  than  the  non- 
specifically  bound  H03~.  These  mechanisms,  influenced  by 

for  the  reduced  mobility  of  soa^”  relative  to  N03". 

Not  treated  here  are  biological  influences  upon  the 
nobility  of  these  anions.  Differential  uptake  of  these 
species  by  plants  affects  the  relative  mobilities  of  these 


anions  as  do  differences  in  the  rates  of  Dicrobially  medi- 
ated injiiobiliaation  and  oxidation. 

Some  models  for  the  chemical  effects  of  acid  deposition 
on  groundwater  and  streamvater  quality  have  included 
estimates  for  these  sinks  and  aouroes.  Selected  chemical 
models  are  discussed  below.  These  models  were  developed  for 
use  in  conjunction  with  hydrologic  models.  Plow  models 
appropriate  to  the  leaching  experiments  in  this  study  are 
discussed  subsequently  to  the  chemical  models. 

Chemical  Models 

Among  the  simpler  models  is  that  of  Reuss  and  Johnson 
(1986) . Among  the  several  cation  exchange  reactions 
occurring  in  soils  that  receive  acid  rain,  these  authors 
consider  only  the  ca^^-Al^*  cation  exchange  reaction  because 
Ca^*  is  quite  similar  to  in  selectivity  with  respect  to 

Al^*  and  these  divalent  species  are  generally  more  abundant 
in  the  soil  than  the  monovalent  bases.  Solution  concen- 
tration of  Al3+  is  controlled  by  acidic  hydrolysis  of 
Al-bearing  minerals.  Equilibria  involving  monomeric  A1 
hydroxides  are  included.  Sulfate  sorption  and  desorption 
influence  the  total  positive  charge  in  solution.  Simple 
carbonate  equilibria  are  also  considered  in  this  model.  For 
specified  inputs  of  8042*,  HO3-,  and  Cl',  taown  initial 
amount  of  Ca^t  adsorbed,  and  partial  pressure  of  CO2, 
equilibrium  concentrations  of  these  various 
calculated. 


The  Model  of  Acidificetion  of  Groundwater  in  Catchments 
(HAGIC)  of  Cosby  et  al.  (1985)  is  conceptually  similar  to 
that  of  Reuss  and  Johnson.  Their  HAGIC  model,  however, 
includes  several  cation  exchange  reactions  and  Al-fluoride 
and  Al-sulfate  monomers.  Sink/source  terms,  due  to  nutrient 
uptake,  microbial  transformations,  or  weathering  release, 
are  included  in  mass  balance  equations  for  the  basic  cations 
and  strong  acid  anions.  Sulfate  sorption  is  presumed 
governed  by  the  Langmuir  equation,  [2]. 

The  Integrated  Lake  Watershed  Acidification  Study 
(ILWAS)  model  was  applied  to  two  watersheds  by  Chang  et  al. 
(1986) . The  chemical  submodel  of  llMhS  includes  the 
inorganic  reactions  1)  cation  exchange  among  WH4*,  Ca^'*',  (f*', 
K*,  Hg2*,  and  Na*,  2)  weathering  of  primary  minerals, 

3)  dissolution  of  Al  and  complexation  with  OH",  F*, 
and  organic  ligands,  4)  simple  carbonate  equilibria,  and 
5)  anion  sorption  and  desorption.  Aluminum  may  either  be  in 
equilibrium  with  gibbsite  or  kinetically  derived  from  it. 
Phosphate  and  S04^~  sorption  and  desorption  are  described 
with  linear  isotherms.  This  model  also  considers  biologic 
processes  including  1)  mineralisation  of  nutrient  species 
via  decomposition  of  organic  matter,  2)  nutrient  uptake, 

3)  respiratory  release  of  CO2,  and  4)  oxidation  of  M. 

Although  the  latter  two  models  provide  considerable 
detail  to  the  inorganic  chemistry,  they  neglect  two  addi- 
tional phenomena  which  may  affect  the  chemistry  of  a soil 
undergoing  artificial  acidification.  Although  Al  is  assumed 


exist  in  solution  as  any  of  several  nonomeric  complexes, 
no  allowance  Is  made  for  adsorption  of  any  form  but  Al^*. 
Implicit  in  the  assumption  that  only  Al^*  is  adsorbed  is 
recognition  of  the  difficulty  in  obtaining  unequivocal 
estimates  of  the  amounts  of  the  different  A1  monomers  on  the 
exchange  complex.  Furthermore,  the  existence  of  solution 


and  sorbed  phase  polymeric  A1  is  ignored.  Bache  and  Sharp 
(1976)  have  demonstrated  that  polymeric  A1  can  be  desorbed, 
and  desorption  increases  with  increasing  concentration  and 
decreasing  pH  of  the  extractant,  These  authors  have  shown 
that  polymerization  reactions  are  important  even  at  a pH  as 
low  as  4.S.  Jardine  et  al.  (1985)  successfully  described  A1 
scrption  on  Ca^'*'  saturated  kaolinite  using  a 2-site  model. 


Aluminum  sorption  occurred  by  equilibrium  cation  exchange 
for  Ca^-t  and  kinecically  controlled  polymerization  on  the 
solid  surface. 

The  effect  of  soil  solution  pH  on  anion  sorption  is 
also  neglected  by  these  models.  If  the  extent  of  sorption 
of  an  input  strong  acid  anion  is  appreciably  Increased  with 
decreasing  pH  of  the  soil  solution,  then  predictive  models 
which  only  consider  the  effect  of  solution  concentration  on 


anion  sorption  would  tend  to  overestimate  the  rate  of  cation 
leaching.  The  chemical  equilibrium  model  used  in  this  study 
deferred  to  the  experimental  complexity  of  rigorously 
treating  the  Al  chemistry.  The  model  does,  however,  attempt 


to  accouht  for  pH  effects  on  anion  sorption  during  soil 


acidification 


sorption  isotherms  derived  from 


axperloents  obtained  from  addition  of  the  corresponding 
strong  acid  to  saaples  of  soil.  Descriptions  of  this  model 
and  of  the  simplified  flow  model  are  presented  in  Chapter  4. 

Transport  Hodels 

For  l-dinenslonal,  steady  flow  through  a parous  medium 
the  time  and  position  dependence  of  concentration  of  a 
reactive  species  is  described  by 


The  several  variables  and  constants  of  equation  [1]  have 
been  previously  defined.  Analytic  solutions  exist  for  a 
variety  of  boundary  conditions  provided  the  term  3S/9t  is 
directly  proportional  to  the  term  3C/3t.  For  a system 
consisting  of  several  competitively  sorbed  solutes  1,  the 
terms  3Si/9t  are  not  directly  proportional  to  the  terns 
3c^/9t.  The  set  of  equations  for  such  a system  are  non- 
linear and  must  be  solved  numerically  using  finite  differ- 
ence or  finite  element  schemes.  Rubin  and  James  (1973)  and 
Valocchi  et  al.  (1981),  for  example,  have  modelled  transport 
involving  multispecies  cation  exchange  by  use  of  the 
Galerkin  finite  element  method. 

schweich  and  Sardin  (1961)  have  presented  an  alterna- 
tive formulation  that  is  computationally  simpler  for  such  a 
system.  Here,  again  considering  l-dimenslonal,  steady  flow, 
the  medium  is  imagined  divided  into  several  equal  length 


Equation  [1] 


segments  or  cells, 
balance  relation  for  each  successive  pair  of  cells  by 
omitting  the  dispersion  term  and  estimating  the  convection 
term  as  a differential.  For  a finite  length  soil  column 
divided  into  J number  of  cells  there  are  J equations  of  the 


where  the  subscript  i refers  to  a particular  cell.  Q is  the 
volumetric  flow  rate,  V is  the  pore  volume,  and  H is  the 
mass  of  soil,  other  symbols  are  as  previously  defined.  By 
these  simplifications  the  partial  differential  equation  [1] 
has  been  reduced  to  an  ordinary  differential  equation. 

During  the  solution  of  these  J equations  numerical  disper- 
sion occurs.  It  is  possible,  therefore,  by  proper  choice  of 
J to  mimic  the  dispersion  phenomenon  that  equation  [13  takes 
into  explicit  account. 

Schultz  and  Reardon  (1983)  simulated  2-dimensional 
transport  in  an  aquifer  using  a mixing-cell  model  combined 
with  an  analytic  solution  for  a tracer  instantaneously 
released  In  a flow  stream.  Their  approach  seems  readily 
adaptable  to  l-dlmensional  flow.  In  this  approximation  the 
flow  medium  was  conceptually  divided  into  an  arbitrary 
number  of  cells.  At  the  beginning  of  a given  time  step  the 
total  solution  mass  of  a solute  within  a cell  was  imagined 


QCj.j  = QC2  + 


« ^1 


concentrated  along  the 


of  the  cell.  During  the 


central  axis 


flow  tins  step  2-dimensional  spreading  occurred  as  the 
solute  initially  concentrated  along  the  axis  moved  several 
cells  downstream.  The  concentration  in  any  cell  was  then 
calculated  as  that  derived  from  its  upstream  parent  plus 
contributions  derived  from  the  two  or  three  longitudinal  and 
transverse  neighbors  of  Che  parent  cell. 

One  of  the  simplest  of  all  published  approaches  is  a 
mixing-cell  model  described  by  Tanji  et  al.  (1967)  for 
predicting  chemical  changes  during  deep  percolation.  In 
this  model  the  profile  is  conceptually  divided  into  an 
arbitrary  number  of  segments.  Kumerical  dispersion  arises 
during  saturated  flow  as  concentrations  along  the  profile 
are  updated  according  to  the  algorithm 

where  ^ is  the  solution  concentration  of  a species  in  the 
segment  during  the  tine  step  of  flow.  The  algor- 
ithm. [6],  is  applied  for  every  species  for  each  increment 
of  flow,  changes  in  solution  concentration  of  the  various 
species,  due  to  chemical  reactions,  are  estimated  using  an 
equilibrium  model  applied  after  each  increment  of  flow. 

In  this  formulation,  more  or  less  dispersion  can  be 
computaticnally  introduced  by  decreasing  or  increasing, 
respectively,  the  number  of  segments  considered.  The 
algorithm.  [6],  assumes  perfect  nixing  between  inflow 
solution,  and  resident  solution.  ^ since  the 


coefficients  of  both  these  terns  are  0.5.  Numerical 
dispersion  can  also  be  effected  by  varying  these  coeffi- 
cients. As  discussed  by  Tanji  (1970J , increasing  or 
decreasing  f in 


a-f)Cj_j_j 


t’3 


effects  decreased  or  increased  numerical  dispersion,  respec- 

Van  onmen  (ISSS)  presented  a theoretical  basis  for 
proper  choice  of  f based  on  experimental  values  of  pore 
water  velocity,  coefficient  of  dispersion,  and  arbitrary 
choice  of  the  length  of  a cell.  The  author's  model,  applied 
to  flow  involving  binary  cation  exchange,  produced  results 
in  reasonable  agreement  with  an  analytical  solution. 


CHAPTER 


MATERIALS  AMD  METHODS 

Colunn  Leaching  E)cparlaentB 
The  soil  used  was  Cecil  series  (Typic  Hapludult]  soil 
taken  from  a forested  site  near  cleoson,  s.c.  (Hydraulic 
conductivity  data  for  Cecii  soil  from  this  site  are  avail- 
able in  cassel,  19S5.)  The  capacity  of  this  soil  to  sorb 
anions  is  well  documented  (Ensminger,  1954;  Kamprath  et  al., 
1956;  Thomas,  1960;  Chang  and  Thomas,  1963).  Columns.  10  cm 
long  by  9.32  cm  diameter,  were  uniformly  packed  with  air- 
dried  topsoii  or  subsoil  materials  to  bullc  densities  of 
1.57  g cn*^  and  1.48  g cn”^,  respectively.  These  columns 
were  then  wet  from  the  bottom  up  with  either  HNO3  or  H2SO4 
solutions  at  pH  3.0,  3.7,  or  4.7.  input  flux  was  maintained 
essentially  constant  and  effluent  fractions  were  collected 
until  approximately  40  pore  volumes  had  been  collected. 

Plow  and  input  concentration  data  for  each  of  the  12  columns 
are  shown  in  Table  3-1. 

concentrations  of  HNO3  corresponding  to  pH  3.0,  3.7, 
and  4.7  solutions  were  calculated  using  the  Davies  eguation 
for  estimating  activity  coefficients.  Standardized  acid 
was  diluted  to  achieve  the  target  concentrations.  The 
diluted  acid  solutions  were  checked  for  final  pH.  The 
concentration  of  each  solution  was  taken  as  the  initial 


appropriate  dilution  factor, 


Preparation  of  K2SO4  aolutions  was  similar,  owing  to  the 
incomplete  second  dissociation  of  the  dibasic  acid,  the 
concentration  In  the  input  solution  was  less  than  the 
normality  of  the  acid.  However,  based  on  a pK^  for  HSO4'  of 
1.96  (Adams,  1971),  H2SO4  was  approximately  97%  dissociated 
at  the  highest  concentration  used.  The  remaining  3%  of  ff* 
tied  up  in  HSOa"  was  presumed  readily  released  once  the  acid 
entered  the  soil.  Hence,  input  concentration  was  taken 
egual  to  the  normality  of  the  acid  solution. 

Samples  of  effluent  were  subsequently  analyzed  for 
1)  Al,  colorimetrically  using  Sriochrome  Cyanine  P (HcLean, 
1965)  f 2)  colorimetrically  by  the  indophenol  blue 

method  (Keeney  and  Nelson,  1982):  3)  ca^t  and  Mg^"’’,  by 
atomic  absorption  spectroscopy;  4)  If*',  by  specific  ion  elec- 
trode; 5)  K*  and  Na'^,  by  atomic  emission  spectroscopy;  and 
6)  5042-,  turbidimetrically.  The  method  of  SO42-  analysis 
was  a modification  of  that  described  by  Rhoades  (1982)  in 
which  reagents  were  proportionally  reduced  to  the  small 
aliquot  of  sample  available. 

Soil  Analyses 

Hater  extractable  and  sorbed  phase  concentrations  of 
the  above  ions  and  Cl~  were  determined  for  the  unaltered 
Cecil  topsoil  and  subsoil  materials.  The  soil  was  prepared 
by  drying,  grinding,  and  sieving  with  a 2 mm  sieve.  Ten- 
gram  samples  of  soil  were  extracted  with  28  ml  H^O  by 


shaking  for  l hr,  centrifuging,  and  vacuum  filtering.  Soil 
samples  including  residual  aqueous  solutions  were  then 
extracted  for  sorbed  ions.  In  each  extraction,  the  soil  and 
residual  solution  were  shaken  with  25  ml  of  the  extracting 
solution  for  1/2  hr,  centrifuged,  and  vacuum  filtered.  This 
extraction  sequence  was  repeated  three  times.  To  remove 
interferences  of  organic  matter  in  the  topsoil  on  304^* 
determinations,  0.1  g of  aotivated  charcoal  was  added  to  the 
samples  prior  to  the  water  extraction.  Calcium,  K*, 
and  Ha'*'  were  extracted  with  IN  NH4OAC  buffered  at  pH  7.0. 
Aluminum  and  NH4*  were  extracted  with  IN  KCl.  chloride  and 
SO42*  were  extracted  with  0.016  M KH2PO4.  Total  acidity  was 
determined  using  the  BaCl2~TEA  method  (Thomas,  1982).  The 
difference  between  total  acidity  and  extracted  Al^'*'  was 
assumed  to  represent  exchangeable  H*.  corrections  were  made 
for  the  retained  aqueous  solutions.  Analysis  for  cl~  was  by 
the  mercuric  thiocyanate  method  (Rhue  and  Kidder,  1983). 

All  other  analyses  were  as  above.  Standards  for  304^' 
(topsoil)  were  similarly  exposed  to  the  activated  charcoal. 
Each  extraction  was  replicated  four  times.  Data  from  these 
extractions  were  used  for  soluble  salts,  total  initial 
amounts  of  each  species  (Table  3-2} , and  for  estimating 
Saines-Thomas  type  selectivity  coefficients  used  in  the 
simulation  model.  This  formulation  for  exchange  selectivity 


coefficient 


[9] 


Where  X is  the  equivalent  fraction  of  a species  or  the 
exchanger,  a is  the  sclution  activity  of  a species,  and  2 is 
the  valence  (Gaines  and  Thomas,  19S3).  The  selectivity 
coefficients  ate  presented  in  Table  4-7,  chapter  4.  The  CEC 
of  the  soil  was  assumed  to  be  equal  to  the  sum  of  extracted 
sorbed  phase  positive  charge,  corrected  for  sorbed  anionic 
charge.  (This  point  is  discussed  in  Chapter  4.) 

Anion  sorption  Isotherms 

Data  for  804^*  sorption  isotherms  with  pH  dependence 
were  obtained  by  equilibrating  10  g samples  of  topsoil  or 
subsoil  with  25  ml  portions  of  H2SO4  of  varied  concentra- 
tions. Soil  samples  and  acid  were  shaXen  for  1 hr,  centri- 
fuged, measured  for  pH,  and  vacuum  filtered.  The  filtrate 
was  then  analyzed  for  804^".  Knowing  the  initial  amount  of 
S04^~  present  (native  plus  added)  and  the  amount  in  solu- 
tion, the  amount  sorbed  was  calculated  by  difference.  Here, 
as  above,  interference  by  organic  matter  in  the  topsoil 
samples  was  eliminated  by  the  addition  of  activated  char- 
coal. Acid  blan)c5  with  charcoal  were  run.  Each  equilibra- 
tion was  duplicated. 

It  was  assumed  Chat  the  concentration  dependence  of 
$04*"  sorption  during  soil  acidification  with  HJSO4  or  HCl 
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could  be  described  with  the  Freundllch  equation  [3].  The 
constants  in  the  Freundiich  egtiation  were  estimated  by  a 
least  squares  procedure  applied  to  the  loqrithnic  fora  of 
the  Freundiich  equation 

In(Si)  - in(ki)  mln(Ci)  [9] 

Similar  isotherms  for  Cl~,  rather  than  N03~,  were 
determined.  The  sorption  behavior  of  non^specifically 
adsorbed  cl"  is  essentially  identical  to  that  of  NO3" 
(Kingston  at  al.,  1972;  Thomas,  1960)  and  Cl"  is  more  easily 
and  reliably  determined. 

Soil  Hinerala  and  Texture 

The  mineralogy  of  the  Cecil  topsoil  and  subsoil  was 
determined  by  x-ray  diffraction  methods  and  the  texture  was 
determined  by  the  pipette  method  (Table  3-2) . Mineralogical 
and  textural  data  were  obtained  from  Liu  (in  preparation) . 

Tracer  Breakthrough  curves 

In  order  to  establish  values  for  input  parameters  of 
the  numerical  model  applied  in  this  study,  tracer  BTCs  were 
needed  for  both  soils.  Chloride  was  used  here  as  the 
tracer.  To  prevent  sorption  or  desorption  of  Cl",  which 
would  otherwise  be  expeoted  for  these  soils,  the  experiments 
were  run  at  constant  pH  and  in  the  presence  of  a preferen- 
tially sorbed  anion.  Air-dried,  uniformly  packed  columns  of 
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DISCUSSION 


weakly  sorbed  anions,  including  NOj")  sorption  by  the  ceell 
topsail  and  subsoil,  respectively.  The  Freundlich  isothems 


Table  4-1.  Freundlich  isotherm  constants  for  3042"  and  Cl" 
sorption  in  Cecil  topsoil  and  subsoil. 
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Data  in  Figures  4-1  and  4-J  show  that  in  both  the 
topsoil  and  subsoil  was  sorbed  to  a greater  eictsnt 

than  was  cl“-  Accordingly,  leaching  of  either  soil  with 
H2SO4  nay  be  expected  to  result  in  a lower  rate  of  cation 
loss  than  leaching  with  KNO3.  As  will  be  presently  shown, 
however,  occurrence  of  anionic  effects  on  cation  leaching  is 
strongly  influenced  by  the  pK  of  the  input  acid. 

Leaching  Experinents 

Application  of  pH  3.0  Acids  to  Colunns  of  Topsoil 

BreaXthrough  curves  for  Ca^^  in  effluent  from  coluans 
of  Cecil  topsoil  that  received  3.0  acid  (Fig.  4-3)  show  a 
slight  retardation  for  the  H2SO4  treatnent  conpared  to  the 
HNO3  treatnent.  Also,  peaX  concentration  was  less  and 
Initial  tailing  was  greater  for  H2SO4  than  for  HNO3.  Total 
amounts  of  each  species  leached  were  estimated  by  inte- 
grating the  BTCs  with  respect  to  numbers  of  pore  volumes. 

Table  4-3  shows  the  ratios  of  the  amount  of  each  basic 
cation  leached  to  the  amoxint  initially  present  in  the  soil. 
These  fractions  were  calculated  from  data  presented  in 
Tables  3-2  and  4-2  and  the  Jtnown  mass  of  soil  per  column. 
Acidic  catione  are  not  included  since  there  were  sources  for 
these  species  during  leaching,  soil  columns  received  inputs 
of  H*  and  Ai^"*"  was  liberated  by  acidic  dissolution  of  soil 
minerals.  Assuming  neither  IC**  nor  Mg^'*'  was  generated  by 
aoidic  dissolution  of  soil  minerals,  the  values  for  the 
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ratios  of  these  apecies  leached  from  the  topsoil  with  respect 
to  the  amounts  initially  present  are  obviously  over- 
estimated. Reference  to  Table  3-2  shows  that  In  the  topsoil 
there  waa  little  of  either  K*  or  Hg^*  initially  present. 
Although  errors  in  measurement  of  extractable  and  Mg2+ 
were  small  in  absolute  value,  such  experimental  errors  (along 
with  experimental  errors  in  measurement  of  effluent  concen- 
tration) probably  account  for  the  apparent  leaching  of 
somewhat  more  K*  and  Hg^*  than  were  initially  present. 

The  ratio  of  the  amount  of  leached  to  that  ini- 

tially extractable  from  the  topsoil  exceeds  one  for  all  acid 
treatments  used.  Mineralisation  of  soil  organic  natter  H nay 
have  been  a source  of  HN4'*’  during  these  leaching  experiments. 
Also,  K*  is  not  a completely  efficient  extractant  for  HH4*  in 
soils  that  exhibit  the  capacity  to  fix  HH4'''  (Keeney  and 
Kelson,  1932).  Although,  in  a study  by  Rich  (1960),  A1 
interlayering  in  vermioulite  greatly  reduced  the  capacity  of 
the  material  to  fix  the  A1  interlayered  vermiculite 

apparently  retained  some  capacity  to  fix  added  NH4'*.  It  is 

capacity  to  fix  either  NHa'*’  or  K*. 

Although  the  initial  quantity  of  Ca^*^  renoved  by  HKO3 
was  greater  than  that  removed  by  H2&O4,  by  40  pore  volumes  of 
effluent,  these  quantities  were  equal.  The  initial  steep 
decreases  in  effluent  ca^***  concentration  were  caused  by 


flushing  soluble  salts  initially  present  in  the  soil  fzom  the 
colunn.  Extensive  tailing  in  the  low  concentration  range 
after  about  10  pore  volumes  suggests  there  was  a slow  release 
of  Ca^'*’.  This  slow  release  could  have  been  due  to  changes  in 
exchange  selectivity  as  the  fraction  of  ca^t  on  the  exchange 
coDplex  approached  zero.  It  could  also  have  been  due  to  slow 
diffusion  in/out  of  regions  of  iramobile  water  within  snail 
aggregates  or  interlanellar  spaces  of  the  aluninunhydroxy 
interlayered  verraiculite. 

Brealcthrcugh  curves  for  Ng^"*'  (Fig.  4-4)  are  sinllar  Co 
those  for  Ca^***.  However,  since  Hg^***  was  less  preferentially 
adsorbed  than  Ca^*  and  was  present  in  a lower  quantity, 
breakthrough  was  earlier  for  Mg^***  ions.  Equal  amounts  of 
Hg^'^  were  leached  by  the  two  treatments. 

Hydrogen  ion  BTCs  (Fig.  4-5)  for  the  HjSO^  treatment 
show  a slightly  slower  rate  of  concentration  Increase  with 
number  of  pore  volumes  than  tor  the  HNO3  treatment.  For  each 
of  the  two  acids,  however,  even  after  40  pore  volumes,  the  H* 
concentration  was  only  about  25%  of  the  input  concentration 
(1.04  jiol(+)/in3) . 

This  buffering  by  the  soil  was  due  not  only  to  cation 
exchange  of  H***  for  adsorbed  Al^'^  or  basic  cation  but  also  to 
K*  reaction  with  a soil  mineral.  Table  4-4  shows  the  sums  of 
leaohed  plus  residually  extractable  Al^*^  for  each  acid 
treatment  compared  to  the  initially  extractable  Al^t.  These 
data  indicate  that  there  was  substantial  A1  production  when 
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leached 


with  pH  3.0 


This  buffering,  therefore,  can  be  accounted  for  on  the 
basis  of  H*  reaction  with  an  Al-bearing  mineral.  If,  for 
example,  kaolinlte  dissolved  and  eguilibriun  was  maintained, 
then  the  ion  product  (Al^*) (8i (OH) 4) (OH") ^ would  be  constant. 
Furthermore,  if  (OH")  were  expressed  in  terms  of  (11+)  and 
the  dissociation  constant  tor  water  and  if  (si(OK)4)  were 
constant,  then  the  activity  ratio  (Al^’’)/ (H*)  ^ would  be 
constant.  Accordingly,  for  leaching  of  the  Cecil  soil  with 
acid,  as  the  output  concentration  of  bases  approached  zero, 
the  effluent  concentration  of  would  have  approached  a 
concentration  less  than  the  input  level.  Effluent  concen- 
trations from  the  Cecil  topsoil  suggest  Chat  dissolution  of 
kaollnice,  in  conjunction  with  the  dissolution  of  quartz 
(controlling  the  activity  of  Si(0H)4),  controlled  the 
activity  ratio  (Al^'*')/ (H+) ^ . If  Al^*  and  were  the  only 
cations  in  the  effluent  and  local  equilibrium  occurred, 
the  output  concentration  of  H***  should  approximately  reach 
0.22  mol(<r‘)/m^  (for  input  pH  3.0  and  based  on  solubility 
constants  adapted  from  Marion  et  al.,  1976),  This  value  is 
very  close  to  observed  values  toward  the  end  (approximately 
40  pore  volumes  of  effluent)  of  Che  leaching  experiments. 

Breakthrough  curves  for  AI^t-  (pig.  4-6)  show  somewhat 
retarded  breakthrough  for  the  H28O4  treatment.  Maximum 
concentrations  exceeded,  however,  the  value  predicted  on  the 
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basis  of  a kaolinlte-quartz  systen.  This  value  is  about 
0.82  mol{+)/n5.)  (Input  concsntratlon,  1.04  moll+j/m^, 
minus  the  above  calculated  final  concentration  of  H*  in  the 
effluent,  0.22  mol (*)/n^ . J These  discrepancies  do  not 
disprove  the  hypothesis  that  Che  kaolinite-quartz  system 
controlled  M solubility.  Although  BTCs  for  A1  are  written 
in  Camis  of  Al^*,  the  quantity  actually  measured  was  total 
Al.  Not  only  did  the  A1  measured  include  Al^t^  also 
several  other  monomers,  polymers,  and  complexes. 

PotassiiuD  breakthroughs  (Fig.  4-7)  show  monotonic 
decreases  in  concentration  throughout  the  runs.  Initial 
rate  for  leaching  loss  was  somewhat  greater  for  the  HNO3 
treatment  but  by  40  pore  volumes  equal  amounts  of  K+  had 
been  leached  by  the  two  treatments.  Breakthrough  curves  for 
the  other  monovalent  cations  are  similar  to  those  for  ir*". 

Sulfate  brea)cthrough  curves  are  shown  in  Fig.  4-s.  The 
BTC  for  SOa*",  HjSOa  treatment,  shows  that  the  input  level 
concentration  of  SOa^"  was  reached  in  the  effluent  attar 
about  15  to  20  pore  volumes.  Attainment  of  SO^^"  sorption 
capacity  coincided  with  only  nominal  increases  in  H’" 
concentration  (Fig.  4-5).  Differences  in  Al^*  concentra- 
tions in  the  effluent  vanished  between  IS  to  20  pore  volumes 
of  effluent  showing  that  sorption  was  no  longer 

effecting  retardation  of  cation  leaching.  The  BTCs  for 
SOa^-  leached  by  HNOj  show  an  early  peak  reflecting  a 
slightly  retarded  loss  of  sOa^"  as  the  watting  front  passed 
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through  the  column.  For  the  HNO3  treatment,  the  rate  of 
304*“  loss  was  low  for  pore  volumes  be/ond  the  early  peat. 
Nitrate  is  a poor  desorbing  agent  for  S042-.  Furthermore, 
SO42-  sorption  was  enhanced  by  the  low  pH  of  the  influent 
solution. 

In  general,  differences  In  cation  BTCs  during  leaching 
produced  by  HjSOe  and  HNO3  at  pH  3.0  can  be  explained  by 
differences  in  the  extent  of  sorption  of  the  input  anion  and 
desorption  of  the  leached  anions.  Retarded  cation  BTCs  for 
the  H2SO4  treatment  imply  that  the  difference,  quantity  of 
SOe*"  sorbed  minus  quantity  of  native  anions  desorbed,  was 
greater  than  the  corresponding  difference  for  the  HNO3 
treatment,  NO3*  sorbed  minus  native  anions  desorbed. 

Effluent  concentrations  of  anions  other  than  SO42*  were  not 
measured , however . 

Application  of  pH  3.7  Acids  to  columns  of  Topsoil 

Somewhat  more  Ca®-'-  was  leached  by  HM03  than  by  H2SO4  in 
the  early  stages  of  acid  application  to  soil  columns 
(Fig.  4-9).  Beyond  about  10  pore  volumes  the  difference 
between  these  treatments  was  negligible.  The  BTCs  for  other 
cations  show  no  difference  between  treatments. 

The  breaXthrough  curve  for  ca^*  leached  by  HKO3  at  pH 
3.0  compared  to  that  for  ca^t  leached  by  the  acid  at  pH  3.7 
shows  that  considerably  more  Ca2+  was  initially  leached  by 
the  lower  pH  treatment  (Fig.  4-10) . Rates  of  loss  beyond 
about  20  pore  volumes,  however,  show  that  more  Ca*"*’  was 


(-*-)  ~ION 


(+)  now 


being  leached  by  the  pH  3.7  treatment.  By  20  pore  volumes 
most  of  the  ca^*^  initially  present  in  the  soil  had  been 
removed  from  the  column  leached  with  pH  3.0  HNO3.  More  Ca^* 
simply  remained  In  the  column  leached  with  pH  3.7  HHO3. 
Comparieon  of  the  BTCs  for  pj{  3,q  3.7,  HNO3 

treatment,  shows  a markedly  decreased  rate  of  loss  at  the 
higher  pH  (Fig.  4*11} . Reduced  leaching  at  this  higher  pH 
occurred  for  all  cations. 

Although  the  effects  of  differences  in  and  H03~ 

sorption  obviously  affected  cation  leaching  rates  by  pH  3.0 
acids,  there  was  apparently  little  effect  of  differences  in 
SO^^"  and  1403“  sorption  on  cation  leaching  rates  by  pH  3.7 
acids.  For  the  Cecil  topsoil,  extremely  low  pH  conditions 
are  required  for  these  anionic  sorption  effects  to  be 
manifested.  Furthermore,  the  retardation  of  cation  leaching 
produced  by  soa^-  anions,  with  respect  to  N03",  was  clearly 
very  much  a secondary  effect  when  compared  to  the  principal 
effect  of  concentration  upon  controlling  cation  leaching 

Application  of  pH  4,7  Acids  to  Columns  of  Topsoil 

Ho  differences  in  cation  leaching  rates  were  observed 
between  treatments,  cation  leaching  rates  were  least  at  pH 
4.7.  Figure  4-10  shows  the  Ca^+  BTC  for  HNO3  input  at  pH 
4.7  compared  to  the  Ca2+  BTCs  for  HNO3  at  pH  3.7  and  3.0. 

An  interesting  behavior  observed  for  pH  4.7  acid  was 
that  the  total  concentration  of  positive  charge  being 
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leached  exceeded,  by  several  times,  the  concentration  of 
input  H*.  considering  the  data  from  Tables  B-1  through  B-7, 
It  can  be  shown  that  the  concentration  of  positive  charge  In 
the  effluent,  late  in  the  leaching  experiments,  was  approxi- 
mately 0,11  iool(  + ) m“^.  The  input  concentration  of  was 
only  about  0.02  mol(+)  however. 

As  earlier  mentioned,  the  nicrobially  mediated  oxida- 
tion of  C,  B,  and  8 in  soil  organic  matter  is  an  internal 
source  of  H*.  Although  the  columns  were  water-saturated, 
this  condition  likely  did  not  limit  potential  microbial 
activity  because  the  input  acid  solutions  were  well- 
oxygenated.  Furthermore,  O2  was  convectively  transported 
through  the  soil  columns.  The  H*  concentration  of  the  acid 
solutions,  however,  may  have  limited  microbial  activity  for 
the  lower  pH  treatments  (Alexander,  1977).  Accordingly,  the 
amount  of  H"*"  internally  produced  by  microbial  oxidation  of 
soil  organic  matter  was  likely  greatest  for  the  pH  4.7 

enhanced  leaching  of  cations  (relative  to  Input  H’'  concen- 
tration) from  columns  that  received  pH  4.7  acids  cannot  here 
be  assessed  independently  of  that  due  to  H*  generated  upon 
net  desorption  of  anions.  A rough  estimate  of  the  relative 
effects  of  these  processes,  however,  may  be  made  from  the 
known  organic  natter  content  of  the  soil  and  assumed  C/N/S 
ratios  and  rates  of  oxidation  of  these  organically  bound 
elements.  The  relation 


Cx  - M(OMV100%)f(X/C)r(4t/t)  (1/EQHT)  (1/V)  CIO] 

where  K is  the  nase  of  the  soil  in  the  coiuan  (g) ; OHt  is 
the  weight  percent  of  organic  natter  present  in  the  soil;  f 
is  the  weight  ratio  of  C to  organic  matter  and  is  assumed  to 
be  1.7;  X/C  is  the  weight  ratio  of  an  organically  bound 
element  to  carbon;  r is  the  annual  rate  of  oxidation  of  the 
element;  iit/t  is  the  fraction  of  a year  over  which  the 
leaching  experiment  was  conducted;  EQWT  is  the  equivalent 
weight  of  the  element  (g  eq“l) ; and  V is  Che  total  volume  of 
acid  solution  that  entered  the  column  (m^)  can  be  used  to 
estimate  the  average  solution  concentration  of  an  anion  (and 
the  corresponding  H*  concentration)  derived  from  the 
microbial  oxidation  of  an  organically  bound  element. 

Assuming  C/N/S  is  I.OO/O.IO/O.OI,  r is  approximately  2%  of 
the  total  per  year  (Alexander,  1977),  and  that  the  leaching 
experiment  lasted  for  about  1/50  of  a year,  then  the  average 
effluent  concentration  of  S04^"  derived  from  oxidation  of 
organic  S in  the  topsoil  would  have  been  about  0.0002  mol(-) 
m"3.  similarly,  the  average  effluent  concentrations  of  M03" 
and  H2CO3  would  have  been  about  0.003  mol(-)  and  0.035 
mol  m“’,  respectively.  Since  H2CO3  is  only  weakly  disso- 
ciated, a concentration  far  less  than  0.035  mol(-)  m“^  HC03“ 
(first  dissociation)  would  have  been  produced  by  oxidation 
of  organic  C.  Considering  the  equilibrium  relation 


(H+) (HC03-)/(H2C03) 


the  concentration  ot  HCO3*  at  pH  S.O  (an  approxinate  value 
for  the  effluent  from  columns  that  received  pK  4.7  acids) 
would  have  been  about  0.003  mol(-)  It  seems  reason- 

able, therefore,  that  the  sum  of  the  average  concentrations 
of  N03~,  and  HC03~  (and  the  corresponding  concen- 

tration) derived  from  microbial  oxidation  of  soil  organic 
matter  was  less  than  about  0.01  molf-}  m~^.  This  concentra* 
tion  is  substantially  less  than  the  difference  between 
concentration  of  effluent  positive  charge  and  input  concen- 
tration of  H*.  0.11  mol(+)  m*5  minus  0.02  mol(+)  m"3.  The 
value  for  effluent  concentration  of  positive  charge  is 
probably  inaccurately  high,  however,  since  all  solution  A1 
was  considered  to  be  Al^*^.  Even  if,  for  the  sake  of 
argument,  the  charge  on  solution  A1  Is  taken  to  be  zero  (an 
unlikely  condition  at  pH  5.0),  internal  production  of  H'*'  by 
oxidation  of  soil  organic  matter  could  not  alone  have 
accounted  tor  the  higher  rate  of  positive  charge  loss  than 
H*  input,  in  this  case  the  solution  concentration  of 
positive  charge  would  have  been  0.04  mol(+)  m'^.  Accord- 
ingly, it  may  be  concluded  that  the  higher  rate  of  positive 
charge  loss  from  columns  of  topsoil  that  rsceived  pH  4.7 
acids  than  rate  of  H^  input  was,  in  large  part,  due  to 
desorption  of  anions. 

Results  for  the  subsoil  leaching  experiments  generally 
paralleled  those  for  the  topsoil.  Breakthrough  curves  when 
pH  3.0  acid  was  applied  show  more  marked  differences  in 


effects  of  the  treatments.  Differences,  though  small,  were 
exhibited  at  pH  3.7, 

Application  of  pH  3.0  Acids  to  Columns  of  Subsoil 

Retardation  of  Ca^*^  leached  by  H2SO4  was  obvious  when 
compared  to  ca^'*’  leached  by  HNO3  (Fig.  «-12).  The  maximum 
Ca^^  concentration  in  column  effluent  was  also  greater  for 
the  HNO3  treatment.  Tailing  in  the  Ca^*  BTCs  occurred  lor 
both  treatments.  Beyond  about  18  pore  volumes  of  effluent 
the  rate  of  Ca^'*'  loss  was  greater  for  the  H2SO4  treatment. 
However,  by  40  pore  volumes  the  quantity  of  Ca^*^  leached  by 
KNO3  still  exceeded  that  leached  by  H2SO4. 

Breakthrough  curves  for  Mg^"*"  (Fig.,  4-13)  are  similar  to 
those  for  Ca^t.  Breakthrough  occurred  later  for  the  H2SO4 
treatment  and  the  concentration  maximum  was  lass.  Tailing 
occurred  and  the  rate  of  Hg^^  leaching  was  greater  for  the 
H2SO4  treatment  beyond  about  14  pore  volumes. 

More  ](*  was  Initially  leached  by  HNO3  than  by  K2SO4 
(Fig.  4-14).  The  BTC  for  K*  leached  by  HNO3  shows  a mono- 
tonic decrease  with  increasing  number  of  pore  volumes.  For 
the  K2SO4  treatment,  however,  the  BTC  shows  a small  peak  at 
about  eight  pore  volumes.  Equal  quantities  of  had  been 
leached  by  the  two  acids  by  40  pore  volumes.  Effluent 
concentrations  of  NH4'^  and  Na*,  for  both  acid  treatments, 
steadily  decreased. 

Beyond  early  minima,  concentrations  of  H'*'  in  the 
effluent  were  greater  for  the  HHO3  treatment  (Fig.  4-15). 
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For  the  two  acids,  H"''  concentrations  in  the  effluent  were  no 
more  than  about  ia%  of  the  input  concentration.  This 
buffering  in  the  subsoil  was  likely  due  to  reaction  with 
gibbsite  controlling  the  activity  ratio,  [Al^"*")  / (fT*")  ^ . if 
Al^*  and  H*  were  the  only  cationic  species  and  the  outflow 
concentration  of  charge  equalled  input  concentration,  the 
effluent  concentration  of  would  be  about  0.14  Boi{»)/n3 
(for  pH  3.0  acid  input  concentration  1.04  Bol(+)/n3  and 
based  on  the  solubility  constant  for  gibbsite  adapted  from 
Marion  et  al.,  1976).  A comparable  value  was  being 
approached  experimentally  by  40  pore  volumes  of  effluent. 

similar  to  the  ion  BTCs,  Al^t  breakthroughs 
(Fig.  4-1$)  show  steady  increases  in  concentration  as 
leaching  progressed.  Outflow  Al^'*'  concentrations  were 
greater  for  the  HNOj  treatment  than  for  HjSO^.  Even  by 
40  pore  volumes  of  effluent,  however,  the  Al^*  concentra- 
tions were  well  below  the  maximum  of  0.90  mol(+)/m5  pre- 
dicted for  the  gibbsite  system. 

The  BTC  for  SO42-  anions  for  the  H2SO4  treatment 
(Pig.  4-17)  shows  substantial  retardation  of  the  input 
SO42-.  Input  concentration  of  SO42-  was  not  reached  in  the 
column  effluent  even  by  40  pore  volumes.  The  304^“  BTC  for 
leaching  with  HHO3  shows  a low,  decreasing  loss  of  S04*'. 
Apparently,  the  low  pH  conditions  effected  by  introducing 
HNO3  at  pH  3.0  into  the  column  largely  offset  the  increased 
tendency  to  desorb  SO42-  anions  by  decreased  concentration. 
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Also,  H03~  is  a poor  competitor  with  504^'  for  soil  sorption 

As  with  the  topsoil,  the  explanation  for  the  compara- 
tive retardation  of  ca2+,  Kg^*,  Al^*,  IT*-,  and  K*  breah- 
throughs  when  the  subsoil  was  leached  with  H2SO4  is  that 
sorption  of  304^*  was  dominant  over  desorptive  losses  of 
other  anions  whereas  for  leaching  with  WD3  no  such  domi- 
nance of  K03“  sorption  occurred  over  3042*  desorption  and 
desorption  of  other  native  anions. 

Application  of  pH  3-7  Acids  to  Columns  of  Subsoil 

For  pH  3.7  acids  differences  in  cation  leaching  rates 
were  observed  between  the  acid  treatments,  Leaching  rates 
of  Ca2*  (Fig.  4-18)  were  initially  greater  for  the  HNO3 
treatment  but  decreased  to  a rate  only  nominally  greater 
than  that  for  the  H2SO4  treatment.  Small  differences  in 
leaching  rates  of  Mg2*  (Pig.  4-19)  were  effected  by  the 
different  acids.  Leaching  rate  of  Hg^*^  remained  somewhat 
greater  throughout  the  run  for  the  HMO3  treatment.  Differ- 
ences between  treatments  for  the  other  cations  were  negli- 
gible. 

Comparison  of  Ca2+  BTCs  for  subsoil  columns  leached 
with  pH  3.0  and  pH  3.7  HKOj  (Pig.  4-20)  shows  substantially 
less  Ca2+  was  leached  by  the  higher  pH  acid  solution.  By  40 
pore  volumes,  however,  effluent  concentrations  of  Ca^*  were 
the  same.  The  effect  of  input  pH  was  more  dramatic  for  Al3+ 
(Fig.  4-21).  Nitric  acid  solution  at  pH  3.7  produced  a 
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steady,  low  concentration  in  the  effluent,  equal  to 
that  from  the  pH  3.0  treatment  up  to  about  12  pore  volumes. 
Beyond  this  point,  Al^"*’  concentration  in  the  effluent  from 
the  column  leached  with  pH  3.0  HN03  substantially  Increased. 
By  40  pore  volumes  of  effluent,  the  Al^*  concentration  from 
the  column  leached  with  pH  3.7  acid  was  only  about  1/20  of 
that  for  the  lower  pH  treatment. 

As  with  the  topsoil,  the  effect  of  SO^^-  sorption, 
relative  to  tJOs"  sorption,  upon  cation  leaching  appears 
limited  to  input  pHs  of  about  3.7  and  lower.  At  such  low  pH 
values  the  retardation  of  cation  leaching  by  interaction 
between  S04^~  anions  and  soil  components  is  small  compared 
to  the  greatly  enhanced  leaching  produced  by  high  concen- 

Appllcation  of  pH  4.7  Acids  to  Columns  of  Subsoil 

No  differences  in  cation  leaching  by  HNO3  and  HjSOa 
solutions  were  observed  for  pH  4.7  solutions,  cation 
leaching  was  least  for  this  pH.  Leaching  patterns  of  ca*"*- 
in  column  effluent  tor  pH  4.7  HNO3  treatment  relative  to  the 
Ca2+  BTCs  for  pH  3.0  and  3.7  acid  solutions  are  given  in 
Fig.  4-20. 

Summing  data  for  positive  charge  concentration  in  the 
effluent  near  the  end  of  the  leaching  experiments  (from 
Tables  B-9  through  B-lS)  gives  an  approximate  positive 
charge  concentration  of  0.08  mol(+)  ra"^.  This  value  exceeds 
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applying  Che  relation  [10]  for  microbial  oxidation  of 
organic  matter  in  Che  subsoil  columns  it  can  be  shown  Chat 
the  average  negative  charge  concentration  in  Che  effluent 
derived  from  oxidation  of  organic  matter  was  about  0.01 
mol(-)  m*3.  As  with  the  topsoil,  desorption  of  anions 
apparently  was  more  important  Chan  organic  matter  decompo- 
sition in  effecting  the  higher  rate  of  loss  of  positive 
charge  than  rate  of  H*  input. 

such  desorption  of  anions  likely  does  not  occur 
naturally  under  forest  conditions,  however.  Atmospheric 
inputs  of  acid  rain  are  commonly  enriched  with  basic  cations 
and  accompanying  anions.  The  acid  solution  reaching 
uppermost  portion  of  the  topsoil  is  further  concentrated 
with  such  ion  species  as  it  passes  through  the  forest  canopy 
and  the  litter  layer.  Under  forest  conditions,  the  solution 
reaching  the  subeoil  is  the  resultant  leachate  from  the 
topsoil  which  has  been  buffered  and  changed  in  ionic 
composition.  In  these  experiments  acid  solutions  wsre 
applied  directly  to  columns  of  topsoil.  Nevsrtheless, 
experimentally  observed  effects  from  the  soil  columns 
demonstrate  the  sensitivity  of  the  soil  to  substantial 
changes  in  the  composition  of  the  soil  solution.  As 
discussed  below,  accelerated  leaching  of  soil  cations  by 
acid  solutions  of  low  pK  may  be  mitigated  by  concurrent  net 
sorption  of  input  anion.  In  this  case,  the  input  of  a 


quantity  of  results  In  the  leaching  loss  of  a less  than 
equivalent  amount  of  native  cations. 

Effluent  concentrations,  by  pore  volume  of  leachate, 
for  all  species  and  all  treatments  are  tabulated  in 
Appendix  B. 

charge  Balance 

A cursory  look  at  data  in  Table  4-2  reveals  that  though 
cation  leaching  increased  as  the  pH  of  the  leaching  solution 
was  decreased,  there  is  no  equivalent  relation  between  input 
quantities  and  leaching  losses  of  soil  cations.  This  was 
especially  true  for  the  subsoil.  For  example,  the  approxi- 
mately 50-fold  increase  in  concentration  of  the  leaching 
solution,  for  pH  1.0  compared  to  4.7  H25O4  solutions, 
resulted  in  lass  than  a 4-fold  increase  in  cations  leached. 
Data  from  Table  4-2  can  be  combined  with  known  amounts  of 
input  acid  and  measured  outflow  amounts  of  504^*  anion  to 
identify  which  treatments  resulted  in  net  retention  or  loss 
of  cations  and  show  the  relation  of  such  retention  or  loss 
to  the  retention  or  loss  of  8042". 

For  H'*'  flow  into  columns  of  air-dried  soil  in  which  no 
reactions  other  than  cation  exchange  and  acid-induced  disso- 
lution of  kaolinite  or  gibbsite  occur,  the  amount  of  poal- 
tive  charge  leached  should  exceed  the  total  amount  of 
input  by  the  amount  of  soluble  sslts  initially  present.  If 
anion  sorption-desorption  occurs,  howsver,  the  difference 
given  as  (positive  charge  out  minus 


soluble  salts)  minus 


70 

positive  charge  input  will  tend  to  be  either  negative  or 
positive  depending  upon  whether  or  not  net  sorption  or 
desorption  of  anions  occurs.  Table  4-5  indicates  those 
treatments  for  which  net  retention  or  loss  of  positive 
charge  occurred.  Values  for  amounts  of  soluble  salts  in  the 
topsail  and  subsoil  were  obtained  from  Table  3-2. 

No  treatment  resulted  in  net  retention  of  positive 
charge  for  oolumns  of  topsoil.  For  the  pK  3.0 
treatment,  however,  the  loss  of  positive  charge  was  small 
compared  to  the  charge  on  A1  leached.  Since  this  value  for 
A1  is  an  ovecestimation,  it  seems  likely  that  positive 
charge  actually  accumulated.  Retention  of  occurred  in 

the  topsoil  for  this  treatment.  Leaching  with  pH  3.7  acids 
caused  minor  net  losses  of  positive  charge.  Leaching  with 
pH  4.7  acids  resulted  in  appreciable  net  losses  of  charge. 
Net  desorption  of  anions  must  have  occurred  since  over- 
estimation  of  charge  borne  by  A1  cannot  account  for  these 
lasses.  Losses  of  SO^^-  were  greatest  at  this  highest  input 
pH. 

Subsoil 

Leaching  Cecil  subsoil  with  pH  3.0  HNOg  or  H2SO4 
produced  a net  retention  of  positive  charge  and  this  effect 
was  especially  pronounced  lor  the  HgSO^  treatment.  Concur- 
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occurred.  The  snail  discrepancy  between  net  retained 
cations  and  net  retained  SOa^'  probably  occurred  because  no 
correction  was  made  for  native  soluble  in  the  estina- 

tion  of  SOi^~  leached.  Snail  net  losses  of  positive  charge 
occurred  at  pH  3.7.  Despite  sone  retention  of  SO42”  for  the 
pH  3.7  H2SO4  treatment,  desorption  and  leaching  of  other 
anions  apparently  occurred  to  a greater  extent.  Leaching 
with  pH  4.7  acids  resulted  in  greater  net  losses  of  positive 
charge  and  net  desorption  of  anions  occurred.  Sulfate 
losses  were  greatest  at  input  pH  4.7. 

From  Table  4-5  it  is  clear  that  the  net  retention  of 
cations  at  low  pH  of  acid  applied  or  net  loss  of  cations  was 
related  to  the  retention  or  loss  of  SOa^'.  As  discussed 
earlier,  S04^~  sorption  occurs  by  specific  adsorption  but 
may  also  occur  by  precipitation  reactions.  The  possible 
occurrence  of  such  precipitation  reactions  during  these  acid 
leaching  experiments  is  examined  below.  Dissolution  of 
baolinite  and  gibbsite  is  also  considered. 

Precipitation  Reactions  of  Sulfate  and  Aluminum  Solubility 

Precipitation  as  a basic  aluminun  sulfate  is  one  of  the 
possible  mechanisms  of  804^'  sorption.  Two  such  compounds 
cited  in  the  literature  as  possibly  contributing  to  804^' 
sorption  are  basaluminite,  Al4(OH)ioS04  ■ SKjO,  and  alunite, 
KAl3(OH) $(804)2.  Investigate  if  either  of  these  com- 
pounds may  have  been  controlling  804^'  solubility  in  these 
leaching  experiments,  effluent  concentration  of  Al3+,  h^, 


K^,  and  S04^~,  corrected  to  activities  using  the  Davies 
equation,  were  combined  to  give  the  ion  products, 
(A13*)^{S0,2)/(h+)10,  and  (K) (Al^*) 3 (SO,^-) (H*) 6 . These 
alternative  lormulations  of  the  ion  products  were  written 
assuming  the  ion  product  for  water,  (H*) (0H“) , to  be  10 
Since  the  measured  value  for  A1  overestimated  Al^*^  by  virtue 
of  the  coexistence  of  several  monomeric  and  polymeric 
species,  values  for  Al^*^  used  in  these  calculations  were 
corrected  to  account  for  a portion  of  the  error  by  con- 
sidering the  mono-  meric  species,  Al^t,  a1(OH)2+.  A1(0B)2+, 
Al(OH)3,  and  AlSOa'',  to  be  in  equilibrium.  Sulfate  con- 
centrations were  corrected  for  the  species,  AlSOa***.  The 
solubility  products  for  basaluminite  and  alunite,  expressed 
in  the  above  forms,  are  10**’^  and  lo”  respectively 
(adapted  from  Adams  and  Rawajfih,  1977J.  The  logarithms  of 
the  calculated  ion  products,  across  a range  of  effluent 
collected,  are  shown  in  Tables  4-6  and  4-7.  calculations 
were  based  on  interpolated  effluent  data  from  four  selected 
columns:  topsoil  that  received  pH  7.0  HNO3;  topsoil  that 

received  pK  3.0  K2^°4''  subsoil  that  received  pH  3.0  HNO3; 
and  subsoil  that  received  pH  3.0  H2SO4. 

Compared  to  literature  values,  the  average  ion  product 
tor  basaluminite  for  each  of  the  selected  columns  was  low, 
indicating  undersaturation  with  respect  to  basaluminite. 
Except  for  the  subsoil  column  leached  with  H25O4,  under- 
saturation  increased  as  leaching  progressed.  Effluent 
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concentrations  froc  this  coluon  gave  relatively  constant  ion 
products  which  were  smaller  by  about  one  order  of  magnitude 
than  the  kgp  for  basaluninite. 

Average  ion  products  for  effluent  from  columns  of 
topsoil  and  subsoil  leached  with  KNO3  were  less  than  the 
literature  value  for  alunite.  On  the  other  hand,  average 
ion  products  for  effluent  from  the  H2SO4  leached  soils 
showed  slight  supersaturation  with  respect  to  alunite. 
Leachate  collected  early  in  the  runs  using  topsoil  shoved 
supersaturation  regardless  of  the  type  of  acid  applied.  The 
degree  of  supersaturation  was  greater  for  the  H2SO4  treat- 
ment. By  40  pore  volumes  effluent  from  topsoil  columns  was 
undersaturated.  An  opposite  effect  occurred  during  leaching 
of  subsoil  columns.  Initial  leachate  for  both  acid  treat- 
ments showed  undersaturation.  As  leaching  proceeded, 
leachate  from  the  column  treated  with  KNO3  became  less 
undersaturated  and  the  leachate  from  the  H2SO4  treated 
column  became  supersaturated  with  respect  to  alunite  after 
only  five  pore  volumes. 

These  calculations  fail  to  determine  conclusively  that 
either  basaluminite  or  alunite  was  controlling  5042-  in 
solution.  Since  it  is  possible  that  chemical  equilibrium 
did  not  exist  locally  in  the  columns  under  the  flow  condi- 
tions that  prevailed  (i.e.,  small  contact  times),  these  ion 
product  calculations  do  not  preclude  the  occurrence  of  these 
minerals  or  their  control  of  SO42-  solubility.  Within  the 


limitations  of  the  estimations  of  Al^*  concentration,  how- 
ever, similar  calculations  for  the  activity  ratio, 
(A1^"'')/(H*)  5,  tend  to  confirm  the  earlier  suppositions  that 
haolinite  and  quart:  controlled  this  ratio  in  the  topsoil 
but  in  the  subsoil  this  ratio  was  controlled  by  gibbsite. 

From  Marion  et  al.  (1976),  assuming  quartz  was  con- 
trolling the  solution  activity,  (HjSO^),  during  acidic 
dissolution  of  kaolinlte,  the  calculated  value  for 
(Al5"*')/(H''')5  is  10^'^®.  Average  values  for  this  ratio  from 
topsoil  leachate  are  for  the  HNO3  treatment,  and 

10^'^®,  for  the  H2SO4  treatment  (Table  4-6).  These  calcu- 
lated values  include  corrections  for  effect  of  ionic 
strength  on  activity  coefficients  and  for  hydrolysis  and  ion 
pair  formation  described  above.  For  both  acid  treatments 
this  ratio  decreased  as  leaching  continued.  High  values 
early  in  the  leaching  process  nay  have  been  due  to  organic 
complexes  or  polymeric  forms,  neither  of  which  were  con- 
sidered in  the  calculations.  Average  values  for  the  subsoil 


H^SO,  treatment  (Table  4-7).  For  gibbsite,  (Al^-*)/ (H*)  ^ is 
10^'®^  (Marion  at  al.,  1976).  By  about  20  pore  volumes  of 
effluent  the  activity  ratio  for  leachate  from  the  HHOj 
treatment  had  closely  approached  this  theoretical  value. 
Values  for  the  H25O4  treatment  decreased  with  pore  volumes 
of  effluent  but  remained  higher  than  10^'®^.  Based  on  these 
calculations,  as  well  as  experimental  BTCs,  it  appears  that 


1-like 


Conservation  of  mass  requires  the  input  solution  from  the 
ith  Qell  to  the  cell  downstream  Co  be  given  by  (13). 

During  the  wetting  phase  the  concentration  in  the  cell 
at  the  wetting  front  was  equal  to  the  input  concentration 
from  its  upstream  neighbor  plus  the  concentration  of  soluble 
salts  dissolved  in  the  cell  as  it  was  saturated.  Exit 
concentrations  during  saturated  flow  were  the  input  concen- 
trations from  Che  terminal  cell. 

Assignment  of  a value  for  f was  made  on  the  Oasis  of  a 
least  squares  comparison  of  simulated  tracer  BTCs  to  experi- 
mental breakthrough  of  Cl*.  When  10  segments  were  used  per 
column,  the  best-fit  value  for  f was  0.S4  for  the  topsoil 
and  0.62  for  the  subsoil.  Figures  4-22  and  4-23  show  tracer 
simulations  compared  to  experimental  cl*  breakthroughs  for 
Cecil  topsoil  and  subsoil,  respectively. 

Chemical  Etmilibriua  Model 

For  each  cell  in  the  flow  domain  and  after  each  flow 
time  step,  the  chemical  model  was  evoked.  The  chemical 
model  required  inputs  of  initial  solution  and  sorbed  phase 
concentrations  of  all  species.  Using  several  i i h-rintii 
relations  and  making  certain  assiujptions  concerning  the 
reactions  involved,  the  model  solved  for  equilibrium  con- 
centrations of  these  species.  The  solution  procedure  is 
described  belcw. 

1.  Knowing  the  initial  solution  and  sorbed  phase 
concentrations  for  all  species,  the  total  quantity,  in 


equivalants,  of  each  epeciee  is  calculated,  for  basic 
cations,  by 


Ti 


fidV  ♦ XiQ 


[14] 


where  is  the  totai  anount  in  equivalents  of  species  i in 
the  cell,  Ci  is  the  molar  concentration,  2i  is  Che  charge  on 
the  species,  ie  the  equivalent  fraction  on  the  exchanger, 
V is  the  volume  of  the  cell,  and  Q is  the  CEC  for  soil  in 
the  cell. 

The  total  quantity  of  an  anion  is  that  in  solution  plus 
that  sorbed.  Neither  total  H'*'  nor  Al^'*'  are  calculated  since 
the  total  quantity  of  each  of  these  changes  during  the 
equilibration  process.  During  acidic  dissolution  of  an  A1 
hearing  mineral,  H*  is  consumed.  Precipitation  of  Al^* 
releases  H***  to  solution.  Anion  sorption  and  desorption,  as 
discussed  below,  respectively  remove  from  or  release  to  the 
soil  solution  one  or  the  other  of  these  acidic  cations.  The 
total  acidity,  TA,  however,  does  not  change  and  it  is  this 
quantity  which  is  calculated.  Total  acidity  is  assumed  to 


TA  . (Ch*  •<-  3C*l3+)V  ♦ + Xai3+)Q  ♦ (Sgo^S-  + SsOj") 

[15] 

The  variables  of  the  first  two  terns  are  analogous  to  those 
previously  defined.  The  variables,  Sso^2-  and  S(jo  are 


the  total  equivalents  of  these  anions  sorbed  in  the  cell. 
Two  points  need  to  be  e>cplalned  here, 

First,  the  model  considers  all  A1  In  the  soil  to  be 
Al^*.  This  is  not  strictly  true  as  previously  discussed, 
but  is  a simplifyinq  assumption  for  the  model.  Second,  it 
is  assumed  that  for  every  equivalent  of  anion  sorbed  there 
is  an  associated  equivalent  of  either  H*  or  Al^*.  For  non- 
speclfically  adsorbed  N03"  this  must  be  true  since  tJOj"  is 
bound  by  Coulombic  attraction  to  a protonated  site.  Fur- 
thermore, for  specifically  adsorbed  so, 2-  either  a proton 
was  consumed  during  positive  charging  of  a site  followed  by 
ligand  exchange  for  H2O  or  was  subject  to  consumption  by 
liberation  of  OH"  aniona  upon  adsorption  at  a neutral  site. 
Desorption  of  SO,^"  anions  requires  displacement  by  OH" 
anions,  presumably  derived  from  HjO,  such  that  H'*'  consumed 
during  adsorption  of  SO, 2-  ig  regenerated  during  the  desorp- 
tion process.  Consequently  any  adsorbed  S0,2"  represents 


the  same  equivalent  quantity  in  potential 
where  SO,2*  adsorption  yields  a negative  ; 
charge  may  be  balanced  by  adscrption  of  ai 
in  a system  which  is  principally  acidic,  t 
such  sites  occupied  by  either  K*  or  Al^*  i 
For  S0,2"  removal  from  the  solution  phase  by  precipitation 
as  a basic  aluminumsulfate  (if  such  occurs)  an  approximately 
equal  quantity  of  Al^*  is  also  precipitated.  If  basalumi- 
nite  is  precipitated,  the  equilibrium  reaction  is 


For  instances 
:e  the  resultant 
cation.  However, 
! fraction  of 
probably  large. 


lOHjO  ■ *l4(OH)ioSO« 


[16] 


precipitated. 


Tbe  otner  basic  aluminumsulfate  cited  is  alunite.  The 


differing  from  equal  quantities  of  SO^*"  and  acid  precipi- 
tated by  the  quantity  of  jr*"  precipitated. 

Parameters  used  in  the  Freundlich  equation  are  given  in 
Table  4-1.  As  applied  in  the  mcdel,  the  values  for  the 
proportionality  constant,  k,  were  scaled  to  allow  for 
expression  of  S in  terms  of  equivalents  per  mass  of  soil  in 
a cell  and  c as  molar  concentration.  A 1/2-lnterval  search 
is  used  to  determine  values  for  sorbed  and  solution  concen- 

cell. 

In  such  a trial  and  error  procedure,  an  initial  esti- 
mate is  made  for  a variable.  This  estimate  lies  halfway 
between  the  established  theoretical  maximum  and  minimum. 
Depending  on  whether  a calculated  value,  using  this  esti- 
mate, is  too  large  or  too  small,  a new  maximum  (equal  to  the 
initial  estimate)  or  a new  minimum  (equal  to  the  initial 


K+  + 3A13-*-  ♦ 


♦ 6H2O  - KAl3(0H)6(S04)2  + 6H+ 


(17] 


variable,  lying  haitway  between  the  upcleted 
made  and  a second  value  is  calculated  from  the  new  estimate. 
This  process  is  repeated,  narrowing  the  possible  range  of 
solutions,  until  satisfactory  agreement  is  obtained. 

3.  solution  concentration  of  positive  charge  is  set 
equal  to  solution  concentration  of  negative  charge  deter- 
mined above.  Knowing  the  concentration  of  positive  charge 
in  solution,  it  is  a fairly  simple  matter  to  determine  the 
estimated  equilibrium  sorbed  and  solution  concentration  for 
the  cations  species.  The  procedure  used  is  to  first  calcu- 
late the  equivalent  fraction  of  ir*"  on  the  exchange  complex 
corresponding  to  a guessed  concentration  of  it*  in  solution. 
These  vslues,  together  with  the  total  quantity  of  each  base, 
then  allow  tor  calculation  of  the  solution  concentration  of 
each  base.  Equilibrium  values  for  solution  concentrations 
are  those  for  which  the  sum  of  equivalents  la  equal  to  the 
total  concentration.  Sorbed  phase  concentrations  are  then 
calculated  by  subtracting  solution  concentrations  from  total 
concentrations  for  each  species.  The  method  is  illustrated 


a.  The  maximum  solution  concentration  of  H***  possible 
is  that  for  which  the  only  cations  present  are  K*  and  Al^*. 
The  minimum  concentration  of  K*  is  zero.  The  search  for 
equilibrium  solution  concentrations  of  the  cations  is  begun 
by  calculating  the  maximum  concentration  of  H*.  This  value 


3Cai3h 


Cmi*  - kCH+3 


[19] 


constant  and  equal  to  that  associated  with  the  mineral 
controlling  solubility.  Concentrations  throughout  this 

that  the  ionic  strength  varies  only  slightly  during  the 
equilibration  process  so  that  activity  coefficients  calcu- 
lated on  the  basis  of  initial  solution  concentrations  in  the 

these  concentrations.  The  substituted  expression,  using 


Ck*  * 3kCHt3  . [21] 

is  solved  using  a Newton-Paphson  numerical  procedure.  The 
concentration  maximum. 


b.  For  this  initial  guess  for  the  concentration  of  H+ 
in  the  solution  phase,  the  corresponding  fraction  of  on 
the  exchange  complex  is  calculated  using 

+ Xai3+  - (TA  - (C„*  + 3Cai3+)V 

- (Sso^2-  + SnOj-))/0  [22] 


X*i3+  - f(XH*)3  [231 

I = Xai3V(Xh+)5  [24] 

this  ratio  having  been  determined  from  extractions  of  the 
unaltered  soil,  If  the  Gaines-Thonas  selectivity  coeffi- 
cient for  exchange  is  assumed  constant  and  the 

solution  activity  ratio,  aM3+/(an+)2,  is  constant,  then  the 
ratio  of  eguivalent  fractions,  f,  must  be  constant.  This 
eguation  is  solved  using  a Nevtcn-Etaphson  procedure.  At 
this  point,  solution  concentrations  of  the  bases  corre- 
sponding to  the  guessed  C(ft-  and  calculated  X)f^  values  are 
determined. 

c.  Eguations  of  the  form 


[25] 


Ci  - Ti/(aiV+((XH*)fi)Q/((YHt)*i/Ti)K(CH+)Si)) 


ai(XH+)*i/(Xi(aH+)*i 


[26] 


is  the  selectivity  coefficient  for  the  specific  exchange 
reaction  and  s are  singie-ion  activity  coefficients  calcu- 
lated from  the  Davies  equation.  Values  for  the  Gaines- 
Thosas  selectivity  coefficients  used  in  the  model  are 
presented  in  Table  4-S.  These  coefficients  were  calculated 
from  the  date  presented  in  Table  3-2. 

In  the  procedure,  corrections  are  made  for  the  ion 
pairs,  caso^®  and  MgS04®,  prior  to  calculating  the  concen- 
tration of  Ca2+  and  Hg2+.  The  method  is  shown  in  the 
computer  program  (Appendix  C) . 


" V " Sa"  * ^=Ca**  * 

* 3Cj^3*  + £27] 

is  made.  Depending  on  whether  the  sum  is  greater  or  less 
than  the  total,  the  estimate  of  C(j+  is  changed  in  accordance 
to  tne  l/2-interv6l  search  scheme  and  steps  b.  and  c.  are 
repeated.  The  search  is  stopped  when  agreement  with  equa- 
tion [27]  is  satisfactory. 

4.  The  equilibrium  solution  concentrations  of  the 
cations  are  used  to  calculate  equivalent  fractions  for  each 
species  cn  the  exchanger  by  equations  of  the  form 


Xi  - (Ti  - 2iCiV)/Q 


(28] 


Fron  initial  values  of  solution  and  sorbed  phase  concentra- 
tions this  conputatlonal  Dodel  will  estimate  equilibrium 
values  for  the  variables. 

Comparison  of  Transport  Modal  Results  to  an  Analytical 
Solution  for  cation  Transport  Involving  Binary,  BooovalenC, 

It  was  of  interest  to  compare  the  performance  of  Che 
mixing-cell  transport  model  to  an  analytical  solution  for  a 
simple  case.  Such  a case  is  binary,  homovalenc,  non- 
preference cation  exchange  occurring  during  steady,  satu- 
rated flow.  Equation  (1),  with  the  term,  P,  set  equal  to 
zero,  describes  such  a flow  system.  In  this  case  3S/9t  is 
directly  proportional  to  3C/3t  and  equation  [1]  reduces  to 


D V |C 

9t  R 3x2  R 3X 


(29] 


The  retardation  factor,  R,  is  given  by 


R - <1  * U ) [30] 

where  Q is  the  cation  exchange  capacity  of  the  soil  and  C( 
is  the  total  concentration  of  cations  in  solution.  Several 
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boundary  conditions.  For  the  particular  syatsn  modelled 
here  (a  column  of  Cecil  topsoil,  as  used  in  this  study)  the 
approximate  boundary  conditions  are  flux  conditions  for  a 
finite  length  column.  The  boundary  and  initial  conditions 
are  given  by 


niS/n  M n 


[31] 
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problem.  If 

A(x,t)  • (C[x,t)-Ci)/(Co-Ci)  [34] 


[3S] 


The  value  for  the  dispersion  coefficient,  D,  was  estimated 
from  a least  squares  fitting  procedure  of  the  analytical 
solution  (for  R ■ l)  to  an  experimental  cl"  Sreahthrough 
curve.  The  length  of  the  soil  column  was  10  cm  and  the 
pore  water  velocity  was  2.42  cm/hr.  The  best-fit  value  for 
D was  1.48  cm^/hr.  Values  for  other  input  parameters  were 
Q - 0.046S7  eq/Xg,  Ct  - 0.10  eq/1,  C = 1.57  kg/1,  and  8 - 

For  the  mixing-cell  transport  model,  the  soil  oolumn 
was  assumed  to  be  divided  into  10  equal  segments.  Ileast 
squares  best-fit  of  simulated  tracer  breakthrough  to  experi- 
mental 01"  breakthrough  was  obtained  for  an  f value  (in  the 
algorithm,  112]  and  [13])  of  0.64.  The  exchange  capacity  of 
each  segment  was  0.003975  eg  and  the  volume  was  0.02264  1. 

For  binary,  homovalent,  non-preference  cation  exchange, 
the  chemical  submodel  was  simple.  Knowing  the  initial 
solution  and  sorbed  phase  concentrations  for  both  cation 
species,  total  equivalents  of  each,  in  a given  segment,  were 
calculated  from  equation  [14],  The  cation  exchange 
equilibrium  expression,  analogous  to  equation  [20],  was 
written  in  terms  of  the  solution  concentration  and  equi- 
valent fraction  of  one  of  the  two  species.  The  resulting 
expression  was  reduced  to  an  equation  containing  only  the 


Cl  = Ct  Ti/(Q  * VCt) 


[36] 


Analogous  variables  have  been  previously  defined.  Knowing 
Ct,  C2  was  found  by  difference.  Furthermore,  Itnowlng  the 
total  equivalents  of  each  species  and  their  solution  con- 
centrations, the  equivalent  fraction  of  each  was  calculated 
from  equations  similar  to  [38]. 

Figure  4-24  shows  elution  curves  for  input  and  leached 
cations  predicted  by  the  analytical  and  mixing-cell  models. 
With  respect  to  the  analytical  solutions,  the  mixing-cell 
model  predicted  excessive  dispersion.  This  behavior  appears 
to  be  characteristic  of  the  model.  A possible  cause  of  this 
discrepancy  may  be  that  the  numerical  model  falls  to  con- 
sider input  and  exit  boundary  conditions  observed  by  the 
analytical  model.  The  locations  of  predicted  breakthroughs, 
however,  were  essentially  the  same  for  both  models. 

Inasmuch  as  the  chemistry  of  the  acid  soil  systems  con- 
sidered in  this  research  probably  exerts  a dominant  influ- 
ence on  ion  leaching  rates,  the  failure  of  the  mixing-cell 
transport  model  to  accurately  describe  dispersive  effects  on 
solute  transport  does  not  limit  the  use  of  this  model  for 
prediction  of  ion  leaching  rates  from  soil. 

Equilibrium  Model  Results 

The  predicted  elution  curve  for  Ca^'*'  leached  from  Cecil 
topsoil  by  pH  3.0  HKQj,  compared  to  the  experimental 
leaching  data,  is  shown  in  Fig.  4-25,  Predicted  peak 
emergence  was  retarded,  occurring  at  about  twice  the  number 
of  pore  volumes  of  effluent  as  experimentally  observed 
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{eight  pore  volumes  compared  to  about  four] . The  height  of 
the  simulated  peak  was  nearly  twice  that  of  the  experimental 
peak.  Also,  the  predicted  decrease  In  Ca^'*'  concentration 
beyond  the  maximum  was  abrupt,  going  to  zero  by  about 
13  pore  volumes.  The  experimental  data  showed  extensive 
tailing. 

The  simulation  model  assumed  all  cation  exchange 
reactions  achieved  equilibrium,  clearly,  this  situation  did 
not  occur.  Konequilibrium  conditions  may  have  been  due  to 
any  of  several  causes  including  1)  variations  in  cation 
exchange  selectivity  with  changes  in  the  sorbed  phase 
composition;  2)  diffusion  in/out  of  stagnant  water  regions 
in  aggregates;  or  3)  cation  exchange  rate  limited  by  the 
kinetics  of  exchange  reaction  or  diffusion  to/from  exchange 

In  simulations  not  discussed  here,  the  author  assumed 
regular  variations  in  exchange  selectivities  as  the  composi- 
tion of  the  exchanger  phase  changed.  The  selectivity 
coefficients  for  K***  replacing  Bi^i'*'  were  arbitrarily  assumed 
to  hyperbollcally  decrease  with  decreasing  equivalent 
fractions  of  on  the  exchanger.  These  simulation 

results,  with  respect  to  those  discussed  in  this  section, 
showed  somewhat  earlier  breakthroughs,  lower  concentration 
peaks  and  curve  tailing.  Although  this  simple  type  of 
exchange  selectivity  variation  did  improve  simulation 
results,  it  did  not  fully  account  for  the  experimentally 
observed  breakthrough  data.  Accordingly,  the  data  may  be 


better  accounted  for  by  a nabU.e/inmobile  water  model  or  a 
two-site,  aquilibrium/lci>’^'<^i'=<  model. 

The  existence  of  isolated  micropore  regions  in  the  soil 
matrix,  through  which  flow  does  not  occur,  tends  to  effect 
an  earlier  emergence  of  concentration  peaks,  lower  maximum 
concentrations,  and  tailing,  Simulations  done  by  Mansell  et 
al.  (19B6)  and  Schulin  et  al.  (1966)  have  demonstrated  such 
effects.  Selim  et  al.  (1976)  used  a two-site,  equilibrium 
and  1st  order  adsorption,  model  to  describe  transport  of  a 
pesticide  in  soil.  Jardlne  et  al.  (1985)  were  able  to 
describe  Al  movement  using  such  a model.  Although  choice  of 
models  is  not  arbitrary,  predicted  results  may  be  similar 
for  each  model.  Hkedi-Kizza  et  al,  (1984),  for  example, 
have  discussed  the  mathematical  equivalence  of  combined 
diffusion/equilibrium  sorption  and  two-site,  equilibrium/ 
kinetic  sorption,  models  for  transport  of  a reactive  solute. 

Figure  4-26  shows  predicted  SO42*  effluent  concentra- 
tions compared  to  experimental  data.  Use  of  the  Freundlich 
sorption  isotherm  led  to  the  prediction  of  a gradual 
decrease  in  804^-  concentration  throughout  the  leaching 
process.  Experimental  data  show  an  early  peak  followed  by 
lower  concentrations  than  predicted.  The  pea)c  in  8042" 
concentration  tended  to  coincide  with  the  minimum  in  IT*’ 
effluent  concentration.  (Compare  to  experimental  data  in 
Pigs.  4-5  and  4-8.)  Conversely,  low  concentration  leaching 
of  SOa^"  beyond  the  peak  occurred  along  with  higher  H* 
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effluent  concentration  values.  The  effect  of  solution  if*' 
concentration  on  anion  sorption  has  been  previously  dis- 
cussed. Apparently,  so^^-  sorption  was  strongly  influenced 
by  solution  concentration  of  H+.  Failing  to  explicitly 
include  this  effect  led  to  underprediction  of  initial  SO42- 
leaching  at  higher  pH  values  and  overprediction  of  SOj^- 
leachlng  at  lower  pH  values,  overprediction  of  8042*,  in 
turn,  contributed  to  overestination  of  solution  concentra- 
tion of  cations. 

This  effect  can  be  seen  in  Fig.  4-27  which  shows 
predicted  and  observed  Al^'*'  concentrations.  Maxinum  Al^* 
concentration  was  predicted  at  about  12  pore  volumes  of 
effluent,  beyond  which  point  the  concentration  slowly 
decreased.  No  similar  behavior  was  experimentally  observed. 
Predicted  breakthrough  of  Al^*  was  retarded  with  respect  to 
experimental  data.  Also,  the  breakthrough  was  abrupt.  The 
retarded  and  steep  predicted  Al^*  breakthrough  occurred 
coincidentally  with  the  predicted  Ca^*  peak  emergence  and 
steep  drop  in  effluent  Ca2+  concentration. 

Figure  4-28  reveals  qualitative  agreement  between 
predicted  and  experimental  leaching  data.  The  occurrence 
and  value  of  the  early  concentration  minimum  were  reasonably 
well  accounted  for  by  the  model.  Predicted  breakthrough 
was,  however,  retarded  and  abrupt,  similar  to  the  predicted 
behavior  for  a13+.  The  apparent  plateau  in  concentration 
beyond  about  12  pore  volumes  of  effluent  exhibits  a slight 
decrease  with  increasing  effluent  volume.  As  with  predicted 
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Al^*  elution,  the  behavior  was  due  to  excessive  prediction 
of  SO^^"  leaching. 

Figure  4-29  shows  predicted  versus  experimentally 
observed  Ca^""'  concentrations  tor  Cecil  topsoil  leached  with 
pH  3.0  H2SO4.  As  was  the  case  for  Cecil  topsoil  leached 
with  pH  3.0  HNO3,  the  predicted  tine  of  peak  emergence  and 
concentration  naxinum  were  much  greater  in  simulation  than 
experiment.  The  model  failed  to  account  for  tailing. 

Despite  these  points  of  discrepancy  between  simulation 
results  and  experimental  data,  the  use  of  anion  sorption 
isotherms  predicted  the  relative  effects  of  NQj'  and  304* 
sorption  on  cation  leaching  rates.  Peak  emergence  was 
predicted  to  occur  later  for  the  topsoil  leached  with  pH  3.0 
H2SO4  than  with  pH  3.0  HHO3.  (Compare  Figs.  4-25  and  4-29.) 

Use  of  pH-Pependent  Anion  Sorption  Isothems 

Contributing  to  the  above  discussed  points  of  dis- 
crepancy between  experimental  data  and  simulation  results 
was  the  inapplicability  of  the  Freundllch  isotherms  for 
describing  anion  sorption  in  a soil  subjected  to  marked 
changes  in  pH.  It  effects  of  solution  concentration  for  H* 
could  be  explicitly  included  in  an  anion  sorption  Isotherm, 
better  model  performance  might  be  expected.  It  is 
proposed  that  a simple  empirical  analog  of  the  Freundllch 
equation  may  offer  substantially  improved  prediction  of 
anion  leaching  behavior.  Thia  analog  is 


Si  - ki(Ci)»(CH*)' 


[37] 
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where  C[^  is  the  solution  concentration  of  n is  a 
constant,  and  the  other  variables  have  been  previously 
defined. 

Kaoprath  et  al.  (1956)  demonstrated  that  a family  of 
sorption  isotherms  may  be  generated  for  S04^~  over  a range 
of  pH  values.  Their  pH-dependent  sorption  isotherms  had 
similar  shapes,  differing  essentially  only  in  increased 
sorption,  at  a given  solution  concentration  of  8042",  with 
increased  acidity  of  the  soil  solution.  It  seemed  reason- 
able to  the  author,  therefore,  that  a single  equation  might 
describe  8042*  (or  any  other  anion)  sorption  over  a range  of 
pH  values  if  the  form  of  the  ecfuation  predicted  increased 
sorption  with  increased  solution  concentration  of  either 
8042*  or  H'*'.  A simple  functional  form  having  this  charac- 
teristic is  equation  (37). 

Experimental  data  used  to  estimate  the  parameters  of 
the  Freundllch  equation  [3]  could  be  fitted  to  equation  [37] 
by  a least  squares  procedure  provided  values  for  K'*'  concen- 
trations were  also  known.  (Anion  sorption  data  is  presented 
in  Appendix  A.)  Thus,  by  such  a fitting  procedure  it  was 
ascertained  that  pH-dependent  8042“  sorption  and  Cl* 
sorption  (considered  representative  of  the  sorption  behavior 
of  all  weakly  sorbed  anions,  including  NO3*)  in  the  topsoil 
could  be  described  by 

®S042-  . 0.29(Cgj,^2.)°-^®(C„+)°-i'‘  [38] 
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Figure  4-30  shows  the  predicted  SO^^-  elution  curve  for 
topsoil  leeched  with  pH  3.0  HNO3  when  equations  [30]  and 
[39]  were  used  in  the  model  compared  to  ejcperinentally 
determined  effluent  concentrations.  Clearly,  explicit 
inclusion  of  the  pH  effect  on  anion  sorption  improved  model 
predictions.  (Compare  to  Fig.  4-26.)  The  early  peak,  due 
to  desorption  at  low  concentrations,  in  effluent  80^2- 
concentration  was  predicted.  Also,  the  low  concentration 
leaching  subsequent  to  this  peak,  due  to  enhanced 
retention  at  higher  H*  concentrations,  was  accurately 
described.  The  position  and  shape  of  the  804^**  peak  were 
not  well-described,  however.  Too  much  desorption  at  low 
concentration  was  predicted.  Better  results  might  be 
obtained  if  better  estimates  of  the  parameters,  k,  m,  and  n, 
were  available.  This  would  entail  inclusion  of  more  data 
sets  (S,  Cl,  Cjft]  in  the  least  squares  procedure.  Further- 
more, it  would  be  desirable  to  have  data  from  sorption 
experiments  in  which  Ci  and  were  varied  in  at  least  more 
than  one  systematic  way. 

Inclusion  of  the  isotherm  equations  [38]  and  [39]  into 
the  equilibrium  model  was  accomplished  by  estimating 
equilibrium  sorption  of  each  anion  at  the  initial  concen- 
tration of  H*  in  a cell  (at  step  2 in  the  model).  Step  3 
estimated  the  equilibrium  solution  and  sorbed  phase 
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concentrations  of  cations,  steps  2 and  3 were  then  repeated 
with  a new  value  for  H*  solution  concentration,  somewhat 
better  estimates  of  equilibrium  sorbed  and  solution  phase 
concentrations  night  have  been  obtained  by  passing  through 
these  steps  several  times  but  at  the  expense  of  more 
computational  time. 

Honecfuilibriun  Model 

Possible  causes  for  the  failure  of  the  equilibrium 
model  to  describe  the  observed  leaching  behavior  have  been 
discussed.  Determination  of  selectivity  variations  for  each 
species  in  a multicomponent  system  is  a formidable  task. 
Determination  of  certain  parameters  which  appear  upon 
introduction  of  mobile/immobile  water  phases  or  kinetically 
controlled  sites  for  cation  exchange  are  often  not  suscep- 
tible to  unequivocal  measurement  (NJcedi-Kizaa  et  al.,  1984). 
Despite  such  difficulties,  the  relative  applicability  of  the 
mobile/immobile  water  and  two-site,  equilibrium/)cinetio, 
models  to  describing  H^-lnducad  cation  leaching  of  the  Cecil 
soil  is  now  considered. 

aggregated  medium  could  be  identically  described  using 
either  a mobile/immobile  water  transport  model  or  a mobile 
water  only  transport  model  that  employed  a dispersion 
coefficient  modified  to  account  for  diffusion  in/out  of 
aggregates  if  the  aggregate  size  or  pore  water  velocity  were 
sufficiently  small.  Therefore,  since  Ci“  ETCs  for  the  Cecil 
topsoil  and  subsoil  could  be  described  by  the  transport 


nodel  witihout  inclusion  of  mobile/ iminobi  le  water  regions 
(see  Figs.  4-22  and  4-22),  the  fitted  parameter  f adequately 
accounted  for  the  contribution  of  diffusion  in/out  of 
stagnant  regions  to  the  overall  dispersion.  Failure  of  the 
equilibrium  model  to  account  for  times  of  peak  concentration 
emergence  or  breakthrough,  however,  indicates  that  some 
fraction  of  the  exchange  complex  was  not  at  instantaneous 
equilibrium  with  the  solution  phase.  7he  cause  of  the 
observed  nonequilibriiui  condition  might  have  been  slow 
cation  exchange  reaction  kinetics  or  diffusion  to/from  sites 
for  exchange.  The  latter  case  may  be  identified  with  the 
occurrence  of  mobile/lmmobile  water  phases.  Here,  however, 
the  volume  of  immobile  water  is  small  and  associated  with 
external  and  Internal  particle  surfaces  for  cation  exchange. 
Accordingly,  diffusion  of  tracer  Cl"  in/out  of  this  volume 
may  have  a small  effect  on  observed  cl~  BTCs.  An  anion 
tracer  is  perhaps  even  largely  excluded  from  this  immobile 
water  volume. 

Noneguilibrium  due  principally  to  reaction  kinetics  or 
diffusion  to/fron  exchange  sites  may  be  accounted  for  by 
evoking  a two-site,  equilibrium/kinetic,  model  for  cation 
exchange.  The  kinetic  portion  of  the  model  includes 


(40) 


where  is  the  equivalent  fraction  of  cation  i on  the 
portion  of  the  exchange  complex  where  exchange  reactione  are 
not  governed  by  instantaneous  equilibrium,  is  the 
solution  concentration  of  i,  is  the  equilibrium  sorption 
coefficient,  and  r^  is  a constant  that  varies  from  cation 
species  to  cation  species, 

Several  new  parameters  appear  in  this  model:  1)  the 

rates  of  reaction  for  kinetically  controlled  cation 
exchange;  and  2)  the  fraction  of  exchange  sites  associated 
with  such  kinetically  controlled  exchange.  fHethods  of 
determining  these  parameters  are  discussed  later  in  this 

Ogwada  and  Sparks  (1986J  have  shown  that  if  the  rate  of 
cation  exchange  is  limited  by  diffusion  to  cation  exchange 
sites,  as  well  as  exchange  reaction  rate,  the  overall  rate 
coefficient,  r,  is  given  by 


r a ^diffusion  ^reaction 
^^diffusion  '*  *^reaction^ 

Using  different  methodologies  for  determining  the  overall 
rate  coefficient,  these  authors  were  able  to  calculate 
values  for  the  diffusion  and  reaction  components  of  the 
overall  rate  coefficient.  According  to  their  data  the 
overall  rate  coefficient  for  K*-Ca2+  exchenge  on  a soil, 
kaollnite,  and  vermioulite,  differed  only  slightly  from  the 
calculated  diffusion  component. 


[40] 
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Therefore,  the  overall  rate  coefficient  ri,  is  approxlmateiy 
determined  by  the  rate  of  diffusion  to/from  exchange  sites 
and  the  equilibrium  sorption  coefficient.  For  diffusion 
in/out  of  interlamellar  spaces  of  veraiculita  should  vary 
appreciably  depending  upon  cation  size.  For  the  aluminum- 
hydroxy  interlayered  material  in  the  Cecil  soil  these 

materials  exhibit  interlamellar  spaces  that  are  propped  open 

cation  species  by  the  methods  discussed  by  ogwada  and  Sparks 
(198S),  for  the  sake  of  illustrating  the  applicability  of 

In  order  to  adapt  equations  of  the  form  [40]  to  model 
application  these  equations  must  be  Integrable.  Within  the 
framework  of  the  mixing-cell  model,  Cj,  in  equation  [40], 


can  be  expressed  in  terms  of  the  total  equivalents,  Ti,  of 
the  species  present  in  a cell  at  any  flow  step,  the  solution 
volume  of  the  cell,  V,  the  CEC  of  the  equilibrium  exchange 
sites,  02'  of  nonequilibrium  exchange  sites, 

Qj'  end  the  equilibrium  sorption  coefficient  for  sites 
governed  by  instantaneous  equilibrium,  giving 

^2  = r^((TiKi^/(V*QiKj^^))-(l*(Q2K.^/V+QlK.^)))X.^)  [44] 

If  and  are  constant  the  equation  can  be  integrated 

X,^  - exp, -Brit,,  [4., 

* “ ^i*'i/(''+QlKi^)  [46] 

8 = 1+(Q2Kj^/(V*0iK^^))  [47] 

and  X,  is  the  initial  equivalent  fraction  of  the 

‘2  initial 

species  adsorbed  onto  type  two,  kinetically  controlled. 

For  the  present  soil  system,  however,  and  are 
derived  from  the  expression  for  the  Gaines-Thomas  selec- 
tivity coefficient  for  H+  exchange  for  cation  i and  are 
given  by 

Kj  = (1/Ki)  (XH’-j/CK+)fi[V('*'H+)*i)  - j-1,2  [48] 


where  Kj  is  the  Gainas-Thomas  selectivity  coefficient. 

These  paraneters  are  not  constant,  of  course.  However, 
since  Xn*j  and  Ck+  nay  be  expected  to  vary  only  slightly 
during  the  adsorption  process  the  integration  can  be  made, 
using  initial  values  for  these  variables,  introducing  only 

Equations  of  the  form  [4S]  nay  be  used  in  a two-site 
model  coupled  to  a mixing-cell  transport  nodel  provided  the 
ri  constants  are  known  and  known.  As  mentioned  above, 

no  attempt  was  made  in  this  study  to  neasure  the  r^  values. 
Neither  was  Qj  estimated,  Instead,  Q2  and  a single  value 
for  ri  were  varied  to  achieve  model  results  in  general 
agreement  with  experimental  results.  Although  such  simula- 
tions serve  to  indicate  the  potential  applicability  of  the 
two-site  nodel  for  explaining  the  experimental  results, 
these  parameters,  all  rfs  and  Qj,  are  susceptible  to 
measurement  (Ogwada  and  Sparks,  I9a«;  Jardine  and  Sparks, 

The  logarithmic  form  of  equation  [45]  can  be  employed 
in  batch  type  experiments  to  determine  the  appropriate 
values  for  the  tis  if  the  CBCs  of  the  equilibrium  and 
kinetic  sites  are  known  and  changes  in  solution  concentra- 
tions of  the  adsorbing  species  with  time  are  measured. 

Since  the  model  assumes  the  existence  of  equilibrium  sites, 
film  diffusion  through  a stagnant  volume  of  water  to  surface 
exchange  sites  must  be  ignored.  The  stirred  batch  method  of 
Ogwada  and  Sparks  (1986)  would  appear  appropriate  for  the 


exchange 


by  interparticular  diffusion  and  reaction  rates. 

aluminohydroxy  vemiculite  (or  interior  exchange  sites  in 
organic  matter)  was  the  causa  of  nonequilibrium,  then  the 

saturating  readily  accessible  exchange  sites  with  a large 
cationic  species  such  as  cetyltrimethylammonlum  bromide 

acidity.  This  adsorbate  was  used  by  Jardine  and  Sparks 
(1984a)  to  block  cation  exchange  at  readily  accessible  sites 
while  leaving  interiamellar  sites  open  for  exchange. 

exchange  in  a soil  that  exhibits  different  types  of  sites 

exchange  is  not  the  same  for  both  types  of  sites.  Jardine 
and  Spar)(s  (1984b),  for  example,  have  demonstrated  that 
K+-Ca2+  exchange  selectivity  for  a Paleudult  was  different 

sites  were  available  to  participate  in  the  exchange 

the  use  of  single  valued  selectivity  coefficients  determined 
using  the  unaltered  soil  (as  previously  described),  quanti- 
fication of  extractable  bases,  Ai,  and  titratable  acidity 
remaining  in  a sail  sample  saturated  with  CTAB  would  allow 


exchange  selectivities  of  the  limited  accessible  exchange 
sites  to  be  estimated.  In  the  subsequent  application  of  the 
tvo-site  model  it  was  assumed  that  exchange  selectivities 
determined  for  the  whole  soil  were  applicable  to  both 
readily  and  limitedly  accessible  exchange  sites. 

Kinetically  controlled  cation  exchange  was  easily 
included  in  the  mixing-cell  transport  and  chemical  equi- 
librium model.  After  each  flow  step,  equilibrium  distribu- 
tions of  adsorbed  and  solution  phases  was  calculated.  Rate 
dependent  adsorption  was  calculated.  Due  to  changes  in 
solution  concentrations  of  cations  brought  about  by  the 
kinetic  exchange  step,  the  equilibrium  solution  and  adsorbed 
phase  concentrations  were  subsequently  updated  to  new 
equilibrium  values.  Figures  4-31  and  4-32  show  predicted 
elution  curves  for  Ca^'*’  and  respectively,  obtained 

from  a two-site  model  compared  to  experimental  data.  The 
model  treated  exchange  at  one-half  of  the  sites  as  kineti- 
cally controlled  (Q2  * 0.00198S).  The  single  rate  coeffi- 
cient used,  ri,  wae  0.2S  hr~^.  The  time  over  which  equation 
t45]  was  integrated  was  (volume  of  water  in  a segment}/ 
(volumetric  flow  cate] . The  model  considered  anion  sorption 
to  be  explicitly  pK-dependent  by  use  of  equations  (36)  and 
(39). 

The  location  and  height  of  the  concentration  peak  for 
Ca^**^  was  adequately  described  using  this  model.  Although 
tailing  was  predicted, the  exact  shape  of  the  concentra- 

tion curve  was  only  qualitatively  described.  Nevertheless, 


(+)  now 


CHAPTER 


CONCLUSIONS 
Leaching  EKperiaients 

Rates  at  which  basic  and  acidic  cations  were  leached 
from  columns  of  Cecil  topsoil  and  subsoil  by  application  of 
simulated  acid  rain  depended  on  the  nature  of  the  acid. 
Losses  of  cations  were  less  if  SO42*,  rather  than  the  weahly 
sorbed  NOj*,  was  the  common  anion.  This  effect  was  more 
pronounced  for  the  Cecil  subsoil,  when  pH  3.0  acid  was 
applied  to  columns  of  subsoil,  net  retention  of  positive 
charge  occurred.  The  small  apparent  net  loss  of  positive 
charge  from  the  topsoil  leached  with  pH  3.0  H2SO4  solution 
was  likely  due  to  overestination  of  charge  borne  by  Al. 
when  pH  3.7  H2SO4  solution  was  applied  to  columns  of  subsoil 
losses  of  Ca^"^  and  Kg2+  were  retarded  with  respect  to  the 
case  when  pH  3.7  HNO3  was  applied.  Ditferencas  in  cation 
leaching  losses  due  to  the  type  of  acid  applied  for  the 
topsoil  were  negligible  at  this  pH.  For  pH  4.7  input  acid, 
however,  there  were  no  differences  in  cation  leaching  losses 
due  to  nature  of  the  acid.  At  this  pH,  however,  there  were 
appreciable  net  losses  of  positive  charge  from  the  topsoil 
and  subsoil  irrespective  of  the  acid  type.  These  losses, 
and  the  smaller  net  losses  at  pH  3.7,  were  largely  attri- 
buted to  desorption  of  native  anions.  It  is  doubtful. 


however,  that  such  behavior  occurs  naturally  in  a forest 
soil  because  the  leaching  solution  is  typically  enriched 

Anionic  effects  on  induced  cation  leaching  froD 
Cecil  soil  appear  to  be  limited  to  pH  values  less  than  3. 
for  the  topsoil,  and  about  3.7,  for  the  subsoil.  At  such 
low  pH  values  anionic  effects  on  cation  leaching  were 
secondary  to  the  principal  effect  of  enhanced  leaching  by 
higher  concentrations  of  ions.  The  exact  range  of  pH 
over  which  these  effects  can  be  seen  in  the  natural  soil 
environment  is  not  known.  Processes  such  as  unsaturated 
flow,  redistribution,  nutrient  cycling,  and  microbial 
activity  could  enhance  or  diminish  these  effects  for  acid 

Leaching  Simulations 

The  eguilibrlun  model  predicted  later  breakthroughs 
than  experimentally  observed.  This  model  also  did  not 
account  for  tailing  of  basic  cation  BTCs.  Such  behavior 
could  have  been  due  to  variations  in  exchange  selectivity 
the  sorbed  phase  composition  of  the  exchanger  changed. 
Also,  diffusion  in/out  of  regions  of  immobile  water  or  non 
instantaneous  rates  of  oation  exchange  could  have  caused 
observed  nonequilibrium, 

The  Freundiich  isotherms  did  not  adequately  account 
pH  effects  on  anion  sorption/desorption.  In  turn,  errors 
estimation  of  anion  concentrations  influenced  predictions 
cation  concentrations.  Explicit  inclusion  of  pH  effects 


anion  sorption  by  use  of  modifies  Freundlich  equations,  [3fl] 
and  [39],  improved  predictions  of  SO42-  concentrations  in 
the  leachate. 

It  was  demonstrated  that  a nonequilibrium  model  for 
cation  exchange  has  the  potential  of  describing  H*-induced 
cation  leaching  from  the  Cecil  soil.  Inasmuch  as  tracer 
BTcs  could  be  described  by  a model  which  did  not  include 
regions  of  immobile  water,  it  was  concluded  that  diffusion 
in/out  of  regions  of  stagnant  water  did  not  significantly 
contribute  to  the  observed  nonequilibrium  conditions. 
Accordingly,  a two-site,  equilibrium/kinetic,  model  for 
cation  exchange  was  applied.  Such  a model  could  account  for 
nonequilibrium  conditions  due  to  either  rate  of  cation 
exchange  or  diffusion  to/from  llmitedly  accessible  exchange 
sites  in  interlamellar  positions  of  the  alumninohydroxy 
interlayered  varmiculite  present  in  the  Cecil  soil. 

By  varying  parameters  for  fraction  of  exchange  sites 
which  exhibit  non-instantaneous  exchange  and  the  rate  of 
cation  exchange,  reasonable  agreement  between  model  end 
experimental  results  was  obtained.  The  model  assumed  equal 
rates  of  cation  exchange  for  all  species.  Better  agreement 
between  model  results  and  experimental  data  could  probably 
be  obtained  by  use  of  different  rate  coefficients  for  each 
cationic  species. 

Although  additional  input  data  would  be  needed  to 
quantitatively  simulate  leeching  from  the  Cecil  soil,  it 
seems  likely  that  reasonably  accurats  predictions  could  be 
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APPENDIX 


COLUMN  LBACHINC  DATA 

The  following  16  Cables  contain  the  complete  column 
leaching  data.  These  tables  fall  into  two  groups.  Tables 
B-1  through  B-8  contain  data  for  the  topsoil  leaching 
experiments  and  Tables  B-9  through  B-16  contain  subsoil 
leaching  data.  Hithin  each  group  the  tables  ace  presented 
in  alphabetical  order  by  name  of  ion  species.  Table  B*l« 
for  example,  contains  effluent  concentration  data  for  Al^ 
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COMPUTER  PROGRAM 


TBe  mixing-call  transport  nodal  and  tha  two-sita 
chanical  subnodel  are  listed  below.  The  chenical  model 
Includes  pH-dependent  anion  sorption.  The  transport  model 
includes  a section  o(  code  to  handle  the  Initial  watting 
phase  (transient  tlow)  of  leaching  for  soil  columns.  Both 


models  needed  for  inclusion  of  nonequilibrium  adsorption  nay 
be  deleted  to  give  code  applicable  to  steady  flow  and 
equilibrium  exchange.  Variables  used  are  defined  for  both 


MIXING-C 


TRANSPORT  MODEL 


C VARIABLE  LIST; 

C ARRAYS- 

C C(I,K)  IS  THE  SOLUTION  CONCENTRATION  ARRAY  WHERE  1-1,9 
C REFERS  TO  SPECIES  (I-l,  NH4;  1-2,  K;  I-I,  HA;  1=4  CA; 

C 1-5,  MG;  I-«,  AL;  1-7,  H;  1-8,  S04;  AND  1-9,  MONOVALENT 
C ANIONS);  AND  K-1,11  REFERS  TO  CELL  NUMBER  ( CELL 
C MO.  1 IS  OUTSIDE  THE  COLUMN  AND  CONTAINS  THE  INPUT 
C SOLUTION) 

C X[I,J,X]  IS  THE  SORBED  PHASE  ARRAY.  DIMENSIONS  I AND  K 
C ARE  AS  ABOVE.  DIMENSION  J-1,2  REFERS  TO  EQUILIBRIUM  OR 
C KINETIC  SITES,  RESPECTIVELY 

C PQWER(I)  STORES  VALENCES  OF  THE  SEVERAL  SPECIES 
C BEGIN(I)  STORES  INITIAL  QUANTITIES  OP  THE  SEVERAL  SPECIES 
C AMOUNT(I)  STORES  RESIDUAL  QUANTITIES  OF  THE  SEVERAL 
C SPECIES  REMAINING  IN  THE  COLUMN  AT  CESSATION  OF  LEACHING 
C SIGMA(I)  STORES  THE  TOTAL  QUANTITIES  OF  IONS  LEACHED 
C CONC  AND  EHTER(I)  ABE  USED  TO  HOLD  RESIDENT 
C CONCENTRATIONS  OF  AND  INPUT  CONCENTRATIONS  TO  A CELL, 

0 RESPECTIVELY,  DURING  A FLOW  STEP 

C CEC(J)  ARE  THE  EXCHANGE  CAPACITIES  FOR  EQUILIBRIUM  AND 
0 KINETIC  SITES 
C SINGLE  VARIABLBS- 

C VOLUME  IS  THE  SOLUTION  VOLUME  OF  A CELL 


SELECTIVITY  COEFFICIENTS  FOB  H EXCHANGE  FOR  NH4 , 
CA,  MG,  AND  AL,  RESPECTIVELY 
R IS  THE  RATIO  XAL/{XH**3) 

COEPl  AND  COEF2  ARE  COEFFICIENTS  OF  THE  SORPTION 
ISOTHERMS  FOR  S04  AND  MONOVALENT  ANIONS 
BXi,  EX2,  EX3,  AND  EX4  ARE  THE  EXPONENTS  IN  THESE 


0PEH(UNIT-10,FILE-'NPTSUS' .STATUS-’ 
OPEN (UmT-12, PILE- 'BTCS08'  .STATUS-' 
READ(10.*)  HSEG,NST0P.CEC(1),CEC(2' 
READ(10.»)  FACT1.FACT2.ZT.  COEFl. 


BBGIN[I)-C(1.2) *VOLU«E»POHER(l)+ 

S X(I,1.2)*(0EC(1)+CEC(2)) 

BEGIN |I) ‘ (FLOAT(NSEG-l) )*BEGIN(1) 

5 CONTINUE 

BEOIN(I)-C(I,2)*VOLUME»POWEB(I)+X(3 

8BGIN(I)=(FLOAr(HSBG-l))*BBGIN(I) 

S CONTINUE 
NCOUNT-0 
MCOUNT-0 


ENTEH(J)-C(J.l) 

1P(K  .Bq!  I)  THEN 
C ( J , K)  -ENTER  (J)  +C  ( J . K) 

CONC-C(J.K) 

C(J. K) -ENTER(J) •FACT1+C0NC*F 
BNTER(J)-ENTEB(J)  *FACT2*CONC 


Irl.L), 


IF(N.  EQ.  100)  NCOUKT^KCOUHtIi 
IF[1JK.  EQ.  100)  MCOUBT=MCOt)NT-H 


EHTSR(J)*C(J,1) 

CONCC(J,K) 

C(J,K)-EHTEH(J)*PACT1+C0>IC«FACT 
ENTER (J) -ENTER (J) *FACT2+C0NC»P 
CONTINUE 
CONTINUE 

PV-FLOAT(I)/FLOAT(NSBG-l) 
PPV-FL0AT(2*I/(HSEG-1) ) 

IF(PPV  .BQ.  2.0*PV)  WRITE(12,100 
S ENTER(1)*1000.0,ENTER(2)*1000.0 
S ENTER(4)*2000.0,ENTER(5)*2000.0 
S BNTER(7)*1000.0,ENTER(8)»2000.0 


CALL  GETRED(Ca.L),X(I,l,t.) 

,L),C(3,L),X(3,1,L),X 

,L).C(5,L).X(5,1.L),X 


AMOUNT(I)=0.0 


) -AMOUNT {I ) *C ( I , J) *VOLOME*POWEB ( I 


COHTIWUE 

WRITE(12!2000)  SIG«A(I) *1000 . 0, AHOUNT(I) 
S BEG1N(1)*1000.0 
CONTINUE 

WBITE(12,*)  NCOUNT 
WRITE(12,*)  MCOUNT 
F0RMAT(X,F5.2,9F7.3) 

FORMAT (X,3F7. 3) 


CHEMICAL  EQUmBHIUM  MODEL 
VAKIA8LE  LIST: 

CNH4,...,  CCL  ARB  THE  SOLUTION  CONCENTRATIONS  OF  THE 
SEVERAL  SPECIES 

CTOT  IS  THE  TOTAL  CONCENTRATION,  IN  EQUIVALENTS,  OF 
CATIONS 

XNH4,...,XK  ARE  THE  EQUIVALENT  FRACTIONS  OF  CATIONS  ADSORBED 
ONTO  EQUILIBRIUM  EXCHANGE  SITES 
XNH42,  . . . ,XH2  ARE  THE  EQUIVALENT  FRACTIONS  OF  CATIONS 
ADSORBED  ONTO  RATE  DEPENDENT  EXCHANGE  SITES 
SORfil  IS  SORBED  S04 
SORB2  IS  SORBED  MONOVALENT  ANIONS 

T1T,...,T5T  ARE  TOTAL  EQUIVALENTS  OF  BASIC  CATIONS 

Tl,...,  T5  ARE  TOTAL  EQUIVALENTS  OP  BASIC  CATIONS 
NH4,.,.,  MS  IN  SOLUTION  AND  ADSORBED  ONTO  EQUILIBRIUM 
EXCHANGE  SITES 

T8  AND  T9  ARE  TOTAL  EQUIVALENTS  OF  S04  AND  MONOVALENT 
ANIONS,  RESPECTIVELY 
TAT  IS  TOTAL  ACIDITY  IN  THE  CELL 

TA  IS  TOTAL  ACIDITY  IN  SOLUTION  AND  ADSORBED  ONTO 
EQUILIBRIUM  SITES 

SKI  AND  SKI  ARE  THE  ION  PAIR  CONSTANTS  FOR  CAS04  AND 
HGS04,  RESPECTIVELY 

CS04E1  AND  CS04B2  ARE  ESTIMATES  OF  THE  CONCENTRATION  OF 
S04  EXCLUSIVE  OF  THE  ABOVE  ION  PAIRS 
A AND  B ARE  COEFFICIENTS  OF  THE  SORPTION  ISOTHERMS  FOR 
S04  AND  MONOVALENT  ANIONS,  RESPECTIVELY 
SIONIC  IS  THE  IONIC  STRENGTH 

GAMKAl,  GAMHA2,  AND  GAMMAS  ARE  ACTIVITY  COEFFICIENTS 
CALCULATED  FROM  THE  DAVIES  EQUATION 
CHECK,  RELATV,  TEST,  EST,  COMP,  AND  DELTA  ARE  USED  TO 
EXAMINE  WHETHER  CERTAIN  CRITERIA  ARE  SATISFIED 
K AND  LL  INDICATE  WHETHER  EQUILIBRIUM  CRITERIA  HAVE  BEEN 
MET  AND  COMMUNICATE  THIS  INFORMATION  TO  THE  MAIN  PROGRAM 
OTHER  VARIABLES  ARE  DEFINED  AND  USED  LOCALLY 


AHED- ( HIGK+ALOH) . 1 


CONTINUE 

HIGH=T9/V 
AKE£>>HIGH/2 . 0 


BELATV=rj/SORB2 
IP(AeS(MlATV)  .L 
IF(REIATV  -LT.  0. 


AHED”(HIGH+ALOW) /2 .0 


AMED=<HI6H+AL0N)/2.0 

CONTINUE 

CTOT=2. 0»CS04+CCL 

SIONIOO.  S*  (CNH4+CK*CNA+4 . 0*CCA+4 . 0*CH6+9. 0*CAL»CH+4 . 0* 
FACTOR— O.S»(SQRT(SIONIC)/(l.O+SQRT(SlONIC))-0.3*SIONIC) 


‘•FACTOR 
•»(4.0*PACTOB) 
•»(9.0*FACTOR) 
L>GA1WA1/GA14HA2 


GAHHAl-1 
GAKHA2-1 
GAHHA3>1 
RATIOl-G 

RATI02-GAHHAI«GAKHA1«GAHHA1/GAHHA3 
RK4“RATI01»R4 
RK5-RATI01*R5 
RK6=RATI02*R6 

DO  WBILE(TBST  -GE.  CTOT/1000.0  .AND.  N .LE.  99) 
CU=(DV*CH-I)/DI 
TEST-ABS(Y) 


AWED=HIGH/2.0 


DO  WHILE<EST  .GE.  ALPHA/100000.0 

XH»(OZ*XH-2)/DZ 

EST=A85(Z) 


CNH4-  (Tl*Rltl»CH/XH)/  (0+BKl*V»CH/XH) 

CK=(T2*RK2*CH/XH)/(5»RK2*V*CH/XH) 

CIIA=(T3*RK3»CH/XH)/(Q+BK3*V*CH/XH) 

CS04E1«CS04/(1.0+GAI01A2»GAMMA2*(SK1*CCA+SK2*CMO) ) 

DO  WHILE(COHP  .GT.  0.001  .AND.  J .LE.  991 
CCA-T4/  (2 .0*V»(1 . 0+GAMKA2»GAHMA2»SKl*CS04El)  tQ*XH*XH/ 
(RK4»CH*CH) ) 

CMG^5/(2.0*V»U-0+GAMMA2»GA«HA2*SK2*CS04E1)»Q*XH*XH/ 

CSO4E2»CSO4/(1.0tGAMMA2»GAMMA2*(SKl*CCA+XK2*CMG) ) 
C0KP»ABS(1.0-CSO4E2/CSO4El) 

CS04E1>CS04E2 

ENDDO 

CA1-RK6*R*CH*CH*CH 

CCA'CCAtSAHHAZ ‘GAHHA2 <SK1*CS04E1*CCA 
CHG-CHG'fGAHHA2  <GAHMA2<SK14CS04E1<CHG 
DELTA-CMH4+CK+CNA+2.0<CMG+3.0<CAL+CH-CT0T 
IF(ABS(DELTA)  .LT.  CTOT/1000.0)  GO  TO  100 
IP(DELTA  .LT.  0.0}  THEM 


AMED- (HIGH+ALOH) /2 . 0 
HIGH-AHED 


AHED- (KIGU* ALOW) /2 . 0 


10  COMTIMUE 
100  CONTIHOE 

XNH4-T1/Q-CNH4 <V/Q 


XNA-T37q-CNA*V/Q 
XCA-T4/Q-2 . 0<CCA<V/Q 
XMG=T5/Q-2.0<CHG<V/Q 


IF(LMN  .EQ.  1)  RETURN 

XNH42-SORBED(TlT,r 

XK2-  S0RBED(T2T,I 
XHA2”  S0RBED(T3T,l 
XCA2=  S0BBED(T4T,I 
XHS2-  S0RBED(T5T,I 
DELTAX-i  .0- (XNH42+XK2-l-XNA2-t 
Z»XN2«'R4XH2<'XH2«XH2-DELTAX 


. DELTAX/100000.0 


EHODO 

XAL2-DELTAX-XH2 

Tl=TlT-Xira42*Q2 

T2=T2T-XK2*Q2 

T3-T3T-XIJA2*Q2 

T4-T4T-XCA2*Q2 

T5-T5T-XMG2»Q2 

TA=TAT-DELTAX*Ca 

200  CONTINUE 


C THIS  FUNCTION  CALCULATES  CATION  ADSORPTION/ 

C DESORPTION  FOB  THE  KINETIC  SITES  ACCORDING  TO 
C EQ[45) 

C 

FUNCTION  S0RBED(T,RK,SH1,SH2, CONCH, CHARG1,CHASG2, VOL, 
$ POW,XSTART,ZZ) 

REAL  T,RK, SHI, SH2, CONCH, CHABG1,CHARG2, VOL, POW,XSTART, 

A=(T*(i!o/RK)*((SH2/CONCH)*»POH))/ 

S (VOL+CHARGl*a.O/RK*((SHl/CONCH)**POH)) 
B=1.0+{CHARG2*(i.O/BK)*((SH2/CONCH)**POW))/ 

S (VOL<-CHABOl*(1.0/RK)  *(  (SHl/CONCH)  **POH)  ) 

C«{1.0/RK)  *(  (SH2/COKCH)**POH) 

SORBED- (1 . 0/B) * (A- (A-B*XSTART) *EXP{-B*2Z/C) ) 

RETURN 


BIOGRAPHICAL  SKETCH 
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